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Abstract 

The matrix integral has many applications in diverse fields. This review article begins 
by presenting detailed key backgroimd knowledge about matrix integral. Then the volumes 
of orthogonal groups and unitary groups are computed, respectively. Applications are also 
presented as well. Specifically, The volume of the set of mixed quantum states is computed 
by using the volume of unitary group. The volume of a metric ball in imitary group is also 
computed as well. 

There are no new results in this article, but only detailed and elementary proofs of existing 
results. The purpose of the article is pedagogical, and to collect in one place many, if not all, 
of the quantum information applications of the volumes of orthogonal and unitary groups. 
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Volumes of orthogonal groups and unitary groups are very useful in physics and mathematics 
||2l|3l. In 1949, Ponting and Potter had already calculated the volume of orthogonal and unitary 
group IITSH . A complete treatment for group manifolds is presented by Marinov [Sl, who extracted 
the volumes of groups by studying curved path integrals lH. There is a general closed formula 
for any compact Lie group in terms of the root lattice IIT^ . Clearly the methods used previously 
are not followed easily. Zyczkowski in his paper ||25]| gives a re-derivation on the volume of 
unitary group, but it is still not accessible. The main goal of this paper is to compute the volumes 
of orthogonal groups and unitary groups in a systematic and elementary way. The obtained 
formulae is more applicable. 

As an application, by re-deriving the volumes of orthogonal groups and unitary groups, the 
authors in ||25H computes the volume of the convex (n^ — 1)-dimensional set D (C”) of the density 
matrices of size n with respect to the Hilbert-Schmidt measure. Recently, the authors in IflTlI give 
the integral representation of the exact volume of a metric ball in unitary group; and present 
diverse applications of volume estimates of metric balls in manifolds in information and coding 
theory. 

Before proceeding, we recall some notions in group theory. Assume that a group G acts on 
the underlying vector space A” via for all ^ G ^ and x G A. Let |x) be any nonzero vector in 
A. The subset G\x) ■= {§^1^) : G is called the Q-orbit of |x) G A. Denote 

--{geG : g|z) = |x)}. 
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We have the following fact: 


g\x)^gix^. 

If x^ = {e}, where e is a unit element of g, then the action of ^ on |x) is called/fee. In this case, 

g\x) ^g/{e} ^g r^g\x). 

Now we review a fast track of the volume of a unitary group. Let us consider the unitary groups 
g =U{n + l),n = 1,2,... To establish their structure, we look at spaces in which the group acts 
transitively, and identify the isotropy subgroup. The unitary group U{n -\-l) acts naturally in 
the complex vector space X = through the vector or "defining" representation; the image 

of any nonzero vector |z) = \xp) G is contained in the maximal sphere of radius 

\\xp\\ since \\xp\\ = || Uxp\\ for U EU{n + \), and in fact it is easy to see that it sweeps the whole 
sphere when U runs though U{n + 1), i.e. the group U{n + \) acts transitively in this sphere. 
The isotropy group of the vector |n + 1) is easily seen to be the unitary group with an entry less, 
that isx^ =lA{n) ®1. Indeed, the following map is onto: 

cp :U{n + l) —> U{n + l)\n + l) = 

If we identity U{n) with ZT(n) 0 1, as a subgroup of ZT(n 0 1), then 

U{n)\n 0 1) = |n 0 1), 

implying that ker cp = U{n), and thus 

U{n + 1)/ ker cp L4{n + l)\n + 1) = = (U{n + 1)/ ker cp) |n 0 1). 

Therefore we have the equivalence relation 

U{n + 1)/U{n) (1.1) 


This indicates that 


That is. 


vol (ZT(n 0 1)) = vol • vol {U{n )), n = 1,2,... 


( 1 . 2 ) 


vol (ZT(n)) = vol X vol (S^) x • • • x vol . 


(1.3) 


We can see this in B- The volume of the sphere of unit radius embedded in R"(n ^ 1), is 
calculated from the Gaussian integral: 
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Now 


/•+°° !' 2 /■+°° 2 
= / / e dtT-dr = / cr„_i(r)e dr, 

JO JS"-i(r) JO 

where e'n-i(?') = fsn-i(r) the volume of sphere S"^^(r) of radius r. Since 

c^n-iir) = cr^-i(l) X r"^\ 


it follows that 

/‘+00 

(V^)" = cr„_i(l) X dr, 

implying that 


/ „ i\ 27T2 


vol(S”-M:= 

V J r{ 2 ) 

(1.4) 

Finally, we get the volume formula of a unitary group: 

„ j n(n+l) 

^ 2^ 7T 2 

v°l ("(")) ^-n rm („-!)!■ 

(1.5) 


2 Volumes of orthogonal groups 

2.1 Preliminary 

The following standard notations will be used IITOl . Scalars will be denoted by lower-case letters, 
vectors and matrices by capital letters. As far as possible variable matrices will be denoted by 
X, Y,... and constant matrices by A, B,_ 

Let A = [fly] be a n X n matrix, then Tr (A) = X])Li«;/ is the trace of A, and det(A) is the 
determinant of A, and ^ over a vector or a matrix will denote its transpose. Let X = [xy] be a 
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m X n matrix of independent real entries Xy's. We denote the matrix of differentials by dX, i.e. 


dX := [dxij] = 

dXffii dXffi2 ' ■ ■ dXffin 

Then [dX] stands for the product of the m x n differential elements 


dxii dxi2 
dX2l dX22 


dxin 

dX2n 


m n 

[dX] (2.1) 

and when X is a real square symmetric matrix, that is, m = n, X = X^, then [dX] is the product 
of the n(n + 1) / 2 differential elements, that is, 

[dX] -fifidxy =nd^y- (2-2) 

7=1 '=7 '>7 

Throughout this paper, stated otherwise, we will make use of conventions that the signs will 
be ignored in the product [dX] of differentials of independent entries. Our notation will be the 
following: Let X = [xij] he a m x n matrix of independent real entries. Then 


[dX] = A-1 dXij 

when [dX] appears with integrals or Jacobians of transformations; 


(2.3) 


[dx] = nr=in"=idx,y 


(2.4) 


when [dX] appears with integrals involving density functions where the functions are nonnegative 
and the absolute value of the Jacobian is automatically taken. 


It is assumed that the reader is familiar with the calculation of Jacobians when a vector of 
scalar variables is transformed to a vector of scalar variables. The result is stated here for the 
sake of completeness. Let the vector of scalar variables X be transformed to Y, where 



Xl 


yi 

X = 


and Y = 



. _ 


. yn. 


by a one-to-one transformation. Let the matrix of partial derivatives be denoted by 


3Y _ rdy, 
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The the determinant of the matrix 


dyi 

dXj 


to Y or Y as a function of X it is written as 
/(Y : X) = det 

and 


( 

'dyf 

V 

dXj 


is known as the Jacobian of the transformation X going 

or [dY] =/(Y:X)[dX], J^O 


J(Y : X) = 


/(X : Y) 


or 1 = /(Y : X)/(X : Y). 


Note that when transforming X to Y the variables can be taken in any order because a permu¬ 
tation brings only a change of sign in the determinant and the magnitude remains the same, that 
is, |/| remains the same where |/| denotes the absolute value of /. When evaluating integrals 
involving functions of matrix arguments one often needs only the absolute value of /. Hence 
in all the statements of this notes the notation [dY] = /[dX] means that the relation is written 
ignoring the sign. 


Proposition 2.1. Let X, Y G R" be of independent real variables and A G be a nonsingular matrix 
of constants. IfY = AX, then 


[dY] = det(Tl)[dX]. 


(2.5) 


Proof. The result follows from the definition itself. Note that when Y = AX, A = [fly] one has 

yi = anxi H-hfl,„x„,f = l,...,n 

where xfs and yj's denote the components of the vectors X and Y, respectively. Thus the partial 
derivative of i/, with respect to Xj is fly, and then the determinant of the Jacobian matrix is 
det(A). □ 


In order to see the results in the more complicated cases we need the concept of a tensor 
product. 


Definition 2.2 (Tensor product). Let A = [fly] G RP^^I and B = [&y] G Then the tensor 

product, denoted by is a pm x qn matrix in formed as follows: 


aiiB 

fll2B ■ ■ 

^IqB 

axiB 

fl22 B ■ ■ 

■ aiqB 

flpi B 

apiB ■ ■ 

UpqB 


( 2 . 6 ) 
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and 


biiA 

fci2A • • 

bin A 

^21A 

^22-^ ’ ■ 

b2nA 


bml-^ ' ' 

bynnA 


(2.7) 


Definition 2.3 (Vector-matrix correspondence). Let X = [x/y] G matrix. Let the j-th column 

of X be denoted by Xy. That is, X = [Xi ,... ,X„], where 


X/ = 



^mj 

Consider an mn-dimensional vector in K"*”, formed by appending Xj,..., X„ and forming a long 
siring. This vector will be denoted by vec(X). That is, 

■ Xi ■ 


vec(X) = 


(2.8) 


X 


n 


From the above definition, we see that the vec mapping is a one-to-one and onto correspon¬ 
dence from to ]R” K'”. We also see that vec(AXB) = (A B^) vec(X) if the product AXB 

exists. 


Proposition 2.4. Let X, Y G be of independent real variables and A G and B G R”^” 

nonsingular matrices of constants. IfY = AXB, then 


[dY] =det(A)Met(B)'"[dX]. 


(2.9) 


Proof. Since Y = AXB, it follows that vec(Y) = (A C) B^) vec(X). Then by using Proposition 12.11 
we have 


/(Y : X) 


det 


aY\ 


I := det 


/ a(vec(Y)) \ 


\3(vec(X)) J 
det(A 0 B^) = det(A 0 !„) det(lm 0 B^) 

det(A)"det(B)'", 


implying that 


[dY] = /(Y : X)[dX] = det(A)” det(B)'"[dX]. 


This completes the proof. 


□ 
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Remark 2.5. Another approach to the proof that [dZ] = det(A)"[dX], where Z = AX, is de¬ 
scribed as follows: we partition Z and X, respectively, as: Z = [Zi,.. .,Z„],X = [Xi,...,X„]. 
Now Z = AX can be rewritten as Zy = AXj for all j. So 


dZ _ ^{Z^,...,Zn) 
dx a(Xi,...,x„) 



A 


( 2 . 10 ) 


implying that 


[dZ] = det 


/ 

’ A 

\ 


A 


V 

A 

/ 


[dX] = det(A)"[dX]. 


Proposition 2.6. Let X,A,B G be lower triangular matrices where A = [fl/y] and B = [bij] are 

constant matrices with ajj > 0, bjj > 0,j = 1,... ,n and X is a matrix of independent real variables. Then 


Thus 


y = X + X^ 
y = AX 

Y = XB 


[dy] = 2'^[dx], 


v;=i 


[d^] = ( 1 [dX]. 


v;=i 


y = AXB ^ [dY]=|n«X 


J 




[dX], 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 


Proof, y = X + X"^ implies that 


xn 

0 ■■ 

■ 0 


Xll 

X21 ■ 

Xftl 


2xn 

X21 ■ 

Xnl 

X21 

X22 • 

■ 0 

+ 

0 

X22 ■ 

■ X2n 

= 

X21 

2X22 ■ 

• Xn2 

Xftl 

X„2 ■ 

Xnn 


0 

0 •• 

Xnn 


Xnl 

Xn2 ■ 

'2.Xnn 


When taking the partial derivatives the n diagonal elements give 2 each and others unities 
and hence [dX] = 2"[dX]. If X = AX, then the matrices of the configurations of the par¬ 
tial derivatives, by taking the elements in the orders (yii,i/2i ■ ■ ■ / J/ni); (3/22/ ■ ■ ■ ,yn 2 )', ■ ■ - lynn and 
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(xii, X 21 ..., Xni); {X22,Xni)} ■ ■ ■} x„„ are the following: 


9 ( 1 / 11 / 3/21 ■ ■ ■ ,ynl) _ 

a(Xii,X 2 l...,X„i) ' d{X 22 ,....X„ 2 ) L I J'---' 

= A[i • • • n - l|i • • • n - 1] = fl„„, 

where A[ii ■ ■ - |;i ■ ■ ■ jy] means that the obtained submatrix from deleting both the d,..., f^rth 

rows and the ji,.. .,jv-th columns of A. Thus 


ay 


d(yu,y2i,--;yni) q 

d(xu,X 2 l,...,Xnl) 

ri 9(y22/--vl/ii2) 

d{X22,-,X„2) 


0 

0 


^ynn 

dXf]f2 - 


A 0 ■■■ 0 

0 A[l|i] ••• 0 


0 0 A[i---n-l|i---n-l] 


We can also take another approach to this proof. In fact, we partition X, Y by columns, respec¬ 
tively, Y = [Yi,..., Y„] and X = [Xi,.. ., X„] . Then Y = AX is equivalent to Yj = AXj, j = 1,..., n. 
Since Y, X, A are lower triangular, it follows that 


yn 

yii 

= A 

Wi 

^21 

/ 

3/22 

= A[i|i] 

X 22 

_ yni _ 


^nl 


ynl 


. ^nl _ 


/ • • • /3/«n — A\l ■ ■ ■ n l|i • ■ ■ n l^Xfm — UfmXfm. 


Now 


[dY] = nidY,] = det(A) det(A[i|i]) • • ■ det(A[i ■ • • n - l|i • ■ ■ n - 1 ]) HldXy] 

1=1 ;=i 




v;=i 
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Next if Y = XB, that is. 


Xn 0 ■ • ■ 

0 


' bn 

0 

■ 0 

X 21 X 22 ■ ■ ■ 

0 


fe2l 

fe22 • 

■ 0 

^nl ^nl 

Xnn 


bnl 

fe,!2 • 

bnn 

xnhi 


0 



0 

3^21 fell + 3^22 fe21 

X22h22 


0 

Xnjbjl 


X 

njbj2 

■ ■ ■ Xyinbnn 


The matrices of the configurations of the partial derivatives, by taking the elements in the orders 
yii; ( 1 / 21 , 1 / 22 ); • • • ; (yni, • • • ,ynn) and Xn; {X 21 , ^ 22 ); ...; {Xni ,..., x„n) are the following: 


Thus 


9yii 

3X11 


= b 


11 , 


9(y3i,y32,y33) 
9(^31, ^32, ^ 33 ) 


^{yiiryii) 

bn fe2i 


fell 0 

^{X21,X22) 

CM 

CM 

0 

_1 


fe21 b22 


T 


' bn 

b2i 

^31 


bn 

0 

0 

0 

b22 

fe32 

= 

bn 

fe22 

0 

0 

0 

^33 


bsi 

bsi 

^33 


9(yni, yn 2 , • • • / y Jin) 
d(^Xnl/ • • • / ■^nn) 


fell 0 

bii ^22 


0 

0 


bnl byi2 • ■ ■ byin 


dY 


3yii 

0 


0 

5(y21,i/22) 


0 0 


0 

0 

^{ynl/Vnl/ ...,ynn') 


Denote by B[ii... .. .jy] the sub-matrix formed by the ii ,... Jji-th rows and ji,... ,;,/-th 

columns of B. Hence 


B[l|l] 0 0 

0 B[12|12] 0 

; 

0 0 
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The whole configuration is a upper triangular matrix with bu appearing n times and ^22 appear¬ 
ing n — 1 times and so on in the diagonal. Also we give another approach to derive the Jacobian 
for Y = XB. Indeed, we partition Y, X by rows, respectively. 



' Vi ■ 


' Xi ■ 

Y = 

L2 

, x = 

X 2 


. _ 


_ X„ _ 


where Yy, Xj are row-vectors. So Y = XB is equivalent to Yj = XjB,j = 1,... ,n. Moreover 

J/ll = = Mlflll/[1/21/1/22] = [^21/^22]B[12|12], ..., 

[i/„i,...,i/„„] = [y„i,...,y„n]B[l...n\l...n]. 


Therefore 


[dY| = nidY,l = det(B[l|l]) det(B[12|12]) ■ ■ ■ det(B[l ■ • ■ n|l ■ ■ ■ nj) nidX;] 
/-I /-I 


YlilX' I ldx|- 


d-i 


We are done. 


□ 


Proposition 2.7. Let X be a lower triangular matrix of independent real variables and A = [fly] and 
B = [by] be lower triangular matrices of constants with ajj > 0, by > 0,/ = 1,..., n. Then 


Y = AX + X^A' 


[dY] = 2" n4 




Y = XB + B^X^ ^ [dY] = 2'M n fe” [dXj. 




(2.15) 

(2.16) 


Proof. Let Z = AX. Then Y = Z -|-Z^. Thus [dY] = 2”[dZ]. Since Z = AX, it follows from 
Proposition 12.61 that 

[dZ] = [dX], 

implying the result. The proof of the second identity goes similarly. □ 


Proposition 2.8. Let X and Y be n x n symmetric matrices of independent real variables and A G 
nonsingular matrix of constants. IfY = AXA"^, then 


[dY] = det(A)"+i[dX]. 


(2.17) 


11 







Proof. Since both X and Y are symmetric matrices and A is nonsingular we can split A and A^ 
as products of elementary matrices and write in the form 


Y = ■ • • E2E^XEJEI ■ ■ ■ 


where Ej,j = 1,2,... are elementary matrices. Write Y = AXA^ as a sequence of transformations 
of the type 

Y^ = E,XEJ,Y2 = E2YiEI...,=^ 

[dYi] = /(Yi : X)[dX], [dY2] = }{Y2 : Yi)[dYi],... 


Now successive substitutions give the final result as long as the Jacobians of the type ]{Y^ : 
Yj._i) are computed. Note that the elementary matrices are formed by multiplying any row (or 
column) of an identity matrix with a scalar, adding a row (column) to another row (column) and 
combinations of these operations. Hence we need to consider only these two basic elementary 
matrices. Let us consider a 3 x 3 case and compute the Jacobians. Let £i be the elementary 
matrix obtained by multiplying the first row by a and E 2 by adding the first row to the second 
row of an identify matrix. That is. 


and 


El 


EiXEj = 


E 2 Y 1 EJ = 


a. 

0 

0 


1 

0 

0 

0 

1 

0 

II 

(N 

1 

1 

0 

0 

0 

1 


0 

0 

1 



UXu 

CiXi3 

1XX21 

X22 

X23 

OCX31 

X32 

^33 


M13 

Ul3 + U 23 

W 33 


Uii Mu + Mi2 

Mu + M21 Mu + M21 + U12 + U22 
U31 M31 + M32 


where Yi = EiXEJ and Y 2 = E 2 Y 1 EJ and the elements of Yi are denoted by Ujfs for convenience. 
The matrix of partial derivatives in the transformation Yi written as a function of X is then 


ax 


0 0 0 0 0 

0 a 0 0 0 0 

0 0 1 0 0 0 

0 0 0 a 0 0 

0 0 0 0 1 0 

0 0 0 0 0 1 
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This is obtained by taking the xij's in the order xn; (x 2 i,X 22 )} (3:3i,a:32, X 33 ) and the Wy's also in 
the same order. Thus the Jacobian is given by 


/(Yi :X) =a^ = = det(£i)3+^. 

Or by definition it follows directly that 

[dYi] = d(a^xii)d(axi2)d(axi3)dx22dx23dx33 = a'^[dX]. 

For an X n matrix it will be Let the elements of Y 2 be denoted by Vij's. Then again taking 
the variables in the order as in the case of Yi written as a function of X the matrix of partial 
derivatives in this transformation is the following: 

’ 1 0 0 0 0 0 ' 

1 1 0 0 0 0 

0 X 2 _ 12 10 0 0 

000100 ’ 

0 0 0 1 1 0 

0 0 0 0 0 1 

The determinant of this matrix is 1 = 1^+^ = det(E 2 )^^^- In general such a transformation gives 
the Jacobian, in absolute value, as 1 = 1”+^. Thus the Jacobian is given by 

/(Y : X) = det(- • • E 2 £i)''+^ = det(A)'’+^. 


We are done. 


□ 


Example 2.9. Let X G be a real symmetric positive definite matrix having a matrix-variate 
gamma distribution with parameters (a, B = > 0). We show that 


det(B)''^ 


1 

T„(a) 


f [dX]det{Xr-^e-^<^^\ 

Jx>o 


^ , n — 1 

Re(a) > — 


(2.18) 


Indeed, since B is symmetric positive definite there exists a nonsingular matrix C such that 
B = CC^. Note that 


Tr (BX) = Tr (C^XC). 


Let 


U = CXC [dli] = det(C)”+^[dX] 
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from Proposition 12.81 and det(X) = det(B) ^det(Lf). The integral on the right reduces to the 
following: 

[ [dX] det(X)'^-'^e-T'-(BX) ^ clet(B)-“ [ [dU] det(!J)“-'^e-T’'(^). 

Jx>0 Ju>0 

But 

[ [dli] det(!J)"-^e-TdLJ) = r„(a) 

Ju>0 

for Re(a) > The result is obtained. 


Proposition 2.10. Let X,Y G skew symmetric matrices of independent real variables and A G 

nonsingular matrix of constants. IfY = AXA^, then 


[dY] =det(A)"-i[dX]. 


(2.19) 


Note that when X is skew symmetric the diagonal elements are zeros and hence there are 
only ”^'2 independent variables in X. 

Proposition 2.11. Let X,A,B G be lower triangular matrices where A = [ap] and B = [by] 

are nonsingular constant matrices with positive diagonal elements, respectively, and X is a matrix of 
independent real variables. Then 

Y = A'^X + X^A ^ [dY] = 2" 

Y = XB^ + BX^ ^ [dY] =2" [dX]. (2.21) 

Proof. Consider Y = A^X + X^A. Premultiply by (A^)~^ and postmultiply by A^^ to get the 
following: 

Y = A^X + X^A ^ (A^)-^YA-^ = (A^)-^X^ + XA“^. 


Let 


Z = XA-i + (XA-I) ^ [dZ] = 2^ [dX] 




by Proposition 12.71 and 


Z = (^A-^yYA-^ [dZ] = det(A)-('’+d[ciy] 
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by Proposition 12.81 Now writing [dY] in terms of [dX] one has 

idri = 2" (pc') idx]=2” (n4) idx] 

since det(A) = 11^1 ^jj because A is lower triangular. Thus the first result follows. The second 
is proved as follows. Clearly, 

y = BX^ + XB^ B-^Y(B^)“1 = B^^X + (^B-^X^ := Z. 

Thus [dZ] = det(B)~(”+^)[dY] and [dZ] = 2” (^II/Li^y/)) [dX]. Therefore expressing [dY] in 
terms of [dX] gives the second result. □ 

Proposition 2.12. Let X G fl symmetric positive definite matrix of independent real variables 

and T = [f/y] a real lower triangular matrix with ty > 0,] = 1,... ,n, and tij,i ^ j independent. Then 

X = T^T [dX] = r [dT], (2.22) 

X = TT [dX] =2" [dr]. (2.23) 

Proof. By considering the matrix of differentials one has 

X = rr^ ^ dx = dr ■ + T ■ dr^. 

Now treat this as a linear transformation in the differentials, that is, dX and dT as variables and 
T a constant. This completes the proof. □ 

Example 2.13. Let X G symmetric positive definite matrix and Re(a) > Show that 

T„(a) := [ [dX]det(X)'*-'^e-'^dx) 

Jx>o 

= 7r^T(a)r(«-i)---r(^-^). (2.24) 


Let T be a real lower triangular matrix with positive diagonal elements. Then the unique 
representation 


X = TT [dX] = T [Yl^l ) [dT]. 


kM 
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Note that 


Tr (X) — Tr {TT~^) — + 1 ^ 2 ) H + (f^i H + ^L) / 

det(X) = det(Tr)=ntJ.. 


When X > 0, we have TT^ > 0 , but tjj > 0,/ = which means that — 00 < tij < 00, i > 

j, 0 < tjj < 00, j = 1 ,..., M. The integral splits into n integrals on tjj's and integrals on ty's, 

i > j. That is. 


But 


T„{a) ^ 



'‘tdtu 


2 




i’dtjj = T 




for Re(a) > = 1,.. .,n and 



Multiplying them together the result follows. Note that 


Re(a) > = h- 


Re(a) > 


n — 1 
2 


2.2 The computation of volumes 


Definition 2.14 (Stiefel manifold). Let A be a m x m{n ^ m) matrix with real entries such that 
A^A = 1,„, that is, the m columns of A are orthonormal vectors. The set of all such matrices A is 
known as the Stiefel manifold, denoted by 0{m,n). That is, for all n x m matrices A, 


0{m,n) = {A G R"'"'" : AM = 1^}, 


(2.25) 


where denotes the set of all n x w real matrices. 


The equation A^A = imposes m{m + l)/2 conditions on the elements of A. Thus the 
number of independent entries in A is mn — m{m + l)/2. 

If m = n, then A is an orthogonal matrix. The set of such orthogonal matrices form a group. 
This group is known as the orthogonal group of m x m matrices. 
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Definition 2.15 (Orthogonal group). Let Bhe an x n matrix with real elements such that B^B = 
In- The set of all B is called an orthogonal group, denoted by 0{n). That is. 


0{n) = {Be : B^B = !„} . 


(2.26) 


Clearly 0{n,n) = 0(n). Note that B^B = 1„ imposes n(n + l)/2 conditions and hence the 
number of independent entries in B is only n^ — n(n + l)/2 = n(n — l)/2. 

Definition 2.16 (A symmetric or a skew symmetric matrix). Let A G If A = A^, then A is 

said to be symmetric and if A'^ = — A, then it is skew symmetric. 

Proposition 2.17. Let V E 0{n) with independent entries and the diagonal entries or the entries in the 
first row ofV all positive. Denote dG = V^dV where V — [vi,... ,v„]. Then 

[dG] = n fl (2.27) 

i=l i=i+l 

= 2”G-i)/2(jgt(l„ + X)-("-i)[dX], (2.28) 

where X is a skew symmetric matrix such that the first row entries of {1„ + X)^^, except the first entry, 
are negative. 


Proof. Let the columns of V be denoted by vi,...,Vn. Since the columns are orthonormal, we 
have {vi,Vj) = 6ij. Then 

{vi,dvj) + {dVi,Vj) = 0 , 

implying {vj,dvj) = 0 since {vj,dvj) is a real scalar. We also have 

{vi,dvj) = -{vpdvi) for i / ;. 


Then V^dV is a skew symmetric matrix. That is. 


dG 


WdL = 


vl 

vl 


[dvi,dv2r- ■ ,dv„] 


{vi,dvi) {vi,dv2) ••• {vi,dvn) 
{V2,dvi) {V 2 ,dV 2 ) ••• {V2,dVn) 


{Vn,dvi) {Vn,dV2) ••• {Vn,dVn) 
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This indicates that 


dG = 


0 {vi,dv2) 

-{vi,dv2) 0 

-{Vi,dVn) -{V2,dVn) 


{Vi,dVn) 

{V2,dVn) 


Then there are only n{n — l)/2 independent entries in G. Then [dG] is the wedge product of the 
entries upper the leading diagonal in the matrix V^dV: 


[dG] 


— a" 

— ^\-=i ^V=!+l 


'vi,dvj) 


This establishes (I2.27I) . For establishing (I2.28I) . take a skew symmetric matrix X, then V = 2{1„ + 
X)~^ — 11„ is orthonormal such that VV^ = Further the matrix of differentials in V is given 
by dV = -2{ln + X)-i ■ dX ■ {!„ + X)-\ i.e. 

dy = -^(i„ + y)-dx-(ii„ + y). 


Thus 


dG = V^dV = -^{1„ + ■ dX • {In + V) 

and the wedge product is obtained 

In + V^y-^ 

V2 ) 

Therefore the desired identity (I2.28II is proved. □ 


[dX] = det (V2(l„ + X)-^j [dXj. 


[dG] = det 


Proposition 2.18. Let X be a n x n symmetric matrix of independent real entries and with distinct and 
nonzero eigenvalues Ai > ■ ■ • > A„ and let D = diag(Ai,..., A„). Let V E 0{n) be a unique such that 
X = VDV^. Then 


where dG = V~^dV. 


[dX] = 


A,--Ail 

) [dD][dG], 


(2.29) 


Proof. Take the differentials in X = 1/D to get 

dx = dv • D • + y ■ dD ■ + y ■ D ■ dy^. 


implying that 


y^ ■ dx • y = y^dy ■ d + dD + d ■ dy^y. 
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Let dY = ■ dX • V for fixed V. Since dG = V^dV, it follows that [dX] = [dY] and 

dY = dG-D + dD + D-dG^ 

= dD + dG-D-D-dG 
= dD+[dG,D] 

where we used the fact that dG^ = —dG and the concept of commutator, defined by [M, N] := 
MN — NM. Clearly the commutator of M and N is skew symmetric. Now 


dyjj = dAj and dyij = {kj — Ai)dgij, i < j = 1,... ,n. 


Then 


n \ / n 


m = yrjdAj.j 

= fn|Ai-A;|j IdDlldG], 

where [dD] = Il/Li dAy and [dG] = n!<; dgij, and the desired conclusion is obtained. □ 


Clearly X = VDV'^, where D = diag(Ai,.. .,A„) (Ai > ■ • • > An),VV~^ = is not a one- 
to-one transformation from X to {D,V) since X determines 2” matrices [±Ui,..., ±y„], where 
Vi,... ,v„ are the columns of V, such that X = VDV^. 


This transformation can be shown to be unique if one entry from each row and column are 
of a specified sign, for example, the diagonal entries are positive. Once this is done we are 
integrating with respect to dG, where dG = V^dV over the full orthogonal group 0{n), the 
result must be divided by 2” to get the result for a unique transformation X = VDV^. 


Remark 2.19. Now we can try to compute the following integral based on Eq. (12.291) : 

/ [dX] = [ ^ I TA.-I I n|A;-A.lfldA,-X / [dG] 

ix>0:Tr(X)=l iAi>-->A„>0 ^ 


A(EA,.-lln|A,-A,indA,x- 


poo ( yi 

n\ Jo 




‘<1 




'O(n) 


Oi(n) 

[dG] 


that is. 


vol 


(DiR'^)) := / [dX] = ^ ^ 

./x>0:Tr(X)=l^ p IlilLtl) ^ f 1 ]) 


, nr 

2 ) 


See below for the notation and vol {0{n)). 
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Proposition 2.20. Let X be a p x n{p ^ n) matrix of rank p and let X = TUJ, where T is a p x p 
lower triangular matrix with distinct nonzero diagonal entries and Ui is a unique n x p semi-orthogonal 
matrix, U^Ui = Ip, all are of independent real entries. Let U 2 be an n x {n — p) semi-orthogonal matrix 
such that Ui augmented with U 2 is a full orthogonal matrix. That is, U = [Jii U 2 ], U^U = LIIU 2 = 
In-p, Lii U 2 = 0. Let Uj be the j-th column of U and duj its differential. Then 


where 


[dX] = 


% 

n-p 

[dr][dUi], 


[dJii] = ^Uj+i (M/vdwy). 


(2.30) 


Proof. Note that 

Ip 0 
0 l„-p 

Take the differentials in X = to get 


WU = 


UJ 


1 

[lil U2] = 

U 2 


UIU2 UIU2 
uLUi um2 


dX = dT-Ul + T- dUf 


Then 


dX-li = dT ■ UlU + T ■ dUlU 

= dT-Ul[Ui, U2] + T-dUl[Ui, U 2 ] 
= [dT + T- dUl -Ui, T ■ dUj ■ U 2 ] 


since UJU^ = Ip, UJL12 = 0. Make the substitutions 


dW = dX • U,dY = dUl ■ Ui,dS = dUj ■U2,dH = T- dS. 


Now we have 


Thus 


dW= [dT + r-dY, dH], 



dfii 

0 ■ 

0 


hi 

0 ■ 

■ 0 


0 

dyi2 • 

■ dyip 

dT + T-dr = 

df2i 

df22 ■ 

0 

+ 

hi 

hi ■ 

. 0 


-dyi2 

0 

■ dy2p 


dtpi 

dfp2 • 

dtpp 


ipl 

ipl ■ 

■ hv - 


—dyip 

-dy2p • 

0 
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Let us consider, for example, the case where p = 2,3 in order for computing the wedge product 
of dT + T ■ dY. Now for p = 2, we have 


dT + T-dY = 

dfii 

0 

+ 

hi 0 

0 

dyi 2 


df2i 

df22 


f2i hi 

_ -dyi2 

0 


dfii fiidyi 2 

dfil — t22dl/l2 df22 + f2idyi2 


Thus the wedge product of dT + T • dY is: 

[dT + T-dY] = dfii A (fiidi/i2) A (df2i — f22dyi2) A (df22 + ^ 21 ^ 1 / 12 ) 
= tiidtii A dyi2 A df2i A df22 = Li [dT] [dY] 

For p = 3, we have 



dfn 

0 

0 


fll 

0 

0 


0 

dyi2 

dyi3 

dT + T-dY = 

df2i 

df22 

0 

+ 

f21 

hi 

0 


-dyi2 

0 

dy23 


df3i 

df32 

df33 


^31 

hi 

^33 


-dyi3 

— dy23 

0 


dfii fiidyi2 fiidyi3 

= df2i — f22dyi2 df22 + f2idyi2 f2idyi3 + f22dy23 

df3i — f32dyi2 — f33dyi3 df32 + f3ldyi2 — f33dy23 df33 + f3idyi3 + f32dy23 

implying the wedge product of dT + T ■ dY is: 

f?if22[dT][dY] = [dr][dY]. 

For the general p, by straight multiplication, and remembering that the variables are only 


dyi 2 / • • • / dyip, dy23/ • • • / dy 2 p, • • •, dy^—^^. 

Thus the wedge product of dT + T ■ dY gives 

ignoring the sign, and 


[dY] = aJL/ {Ui,duj). 

Now consider dH = T ■ dS. Since dS is a p x {n — p) matrix, we have 

[dH] =det(T)"-P[dS] ^ [dS], 
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The wedge product in dS is the following: 


[dS] = A’Lp+; {Ui,duj). 


Hence from the above equations. 


[dX] = [dW] = A[dT + TdY] A [dH] 



Now 


[dy][dS] = aJL/ {Ui,duj) aJLi A^^p^j{ui,duj). 


Substituting back one has 


[dX] = \tjjrn[dT]A^^,Alj^,{ui,dUj) 



which establishes the result. 


□ 


If the triangular matrix T is restricted to the one with positive diagonal entries, that is, tjj > 
0,] = 1,... ,p, then while integrating over T using Proposition l2.20[ the result must be multiplied 
by 2P. Without the factor 2P, the tjj's must be integrated over —oo < tjj < oo,j = l,...,p. If 
the expression to be integrated contains both T and U, then restrict tjj > 0 ,; = 1 ,..., p and 
integrate U over the full Stiefel manifold. If the rows of U are mi, ..., Up, then ±Mi, ..., ±Mp give 
2P choices. Similarly tjj > 0,tjj < 0 give 2P choices. But there are not 2^P choices in X = TU. 
There are only 2P choices. Hence either integrate out the tjj's over —oo < tjj < oo and a unique U 
or over 0 < tjj < oo and the U over the full Stiefel manifold. 

Proposition 2.21. Let Xi be an n x p{n ^ p) matrix of rank p of independent real entries and let 
Xi = UiTi where Ti is a real p x p upper triangular matrix with distinct nonzero diagonal entries and 
Ui is a unique real n x p semi-orthogonal matrix, that is, UJUi = Ip. Let U = [Hi U 2 ] such that 
U^U = In, U 2 U 2 = In-p/ L[JL [2 = 0. Let Uj be the j-th column ofU and duj its differential. Then 



(2.31) 


where 


P n 


[diii]=n n 


M '=/+i 


22 





Proposition 2.22. If Xi,Ti and Ui are as defined in Proposition 12.211 then the surface area of the full 
Stiefel manifold 0{p,n) or the total integral of the wedge product 0{p,n) is 

given by 



where 

Tp{a) = Ti^T{a)T - 0 • • ■ 

for Re(a) > 

Proof. Note that since Xi is n x p, the sum of squares of the np variables in Xi is given by 

Tr(X|Xa) = f;f 4 

Then 



Tr(xTXi) 




_ yn yV y2 
5 L^i=l L^j=\ 


n V 

ri n 


i=i,=i 



HE 
71 2 


by direct evaluation of the exponential integrals. Make the transformation as in Proposition l2.21l 

Xi = UiTi ^ X[Xi = Ti^Ti, 


where T-[ is a real p x p upper triangular matrix with distinct nonzero diagonal entries and Ui is 
a unique real n x p semi-orthogonal matrix—iJi = Ip, implying 

Tr(X[Xi)-Tr(T/Ti) = X:f|. 

N; 


Note that dXi is available from Proposition 12.211 Now 
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But for 0 < tjj < CO, —00 < ty < oo(/ < j) and Ui unrestricted. 



;g-D,;f?;[dri] =2-prpg) 


observing that for ; = 1 ,..., p, the p integrals 



,n~j 


'iidtjj = 2- 


‘ri^ 



,n> i -1, 


and each of the p(p — l)/2 integrals 



= \fn, X < ]. 


Thus the result that follows. 


□ 


Theorem 2.23. Let X be a full-ranked and n x n matrix of independent real entries and let X = UT, 
where T is a real n x n upper triangular matrix with distinct nonzero diagonal entries and U is a unique 
real orthogonal matrix. Let Uj be the j-th column of U and duj its differential. Then the volume content of 
the full orthogonal group 0{n) is given by 



Proposition 2.24. Let Xbe a p x n{p ^ n) matrix of rank p and let X = TUJ, where T is a p x p lower 
triangular matrix with distinct positive diagonal entries tjj > 0,/ = 1 ,..., p and Ui is a unique n x p 
semi-orthogonal matrix, UjUi = Ip, all are of independent real entries. Let A = XX"^ = TT"^. Then 

[dX] = 2^Pdet(A)2-^[dA] A”^j^^(uj,duj}. 


Proof. Since A = TT^, it follows fhaf 


IdA] = I” I I [dT], 




i.e. 


Buf 


[dT]=2-p{YlC^] [dA], 




[dX]= ni^;. 


\n-] 


[dr][dLli], 
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where 


[dill ]=n n 

M ;=/+i 

Note that det(A) = det(T)^ = n]Li The desired conclusion is obtained. □ 


Proposition 2.25. Let X be a m x n matrix of rank m{m ^ n), T = [tyj.] a mxm lower triangular 
matrix with tjj > 0,j = l,...,m and L a n x m matrix satisfying L^L = 1^, where the matrices are of 
independent real entries. Then show that, if X = TL^, then 


[dX] = [dT][dL] 


(2.34) 


where 

m n 

dT=n n 

j=i i=j+i 

Ij is the j-th column ofL= [L Li] G 0(n); dZ, the differential of the i-th column of L. 


Proposition 2.26 (Polar decomposition). Let X be a m x n{m ^ n) matrix, S a m x m symmetric 
positive definite matrix and L a n x m matrix with L^L = all are of independent real entries. Then 
show that, if X = VSL'^, then 

/1 \ - -1 

[dX] - (^-J (det(S))"^ [dS][dL] (2.35) 

where dL is defined in Proposition 12.251 


Proof. Now if X = y/SL^ and L^L = 1,„, then XX^ = S. By Proposition I2.25[ we have X = TTP, 
where T = [tjj^] is am x m lower triangular matrix with tjj > 0,] = 1,... ,m and L an x m matrix 
satisfying L^L = 1^. Denote L = [L Li] G 0{n). Hence 

[dx]=(frr) 


It also holds that S = TT^ implies 


[ds]=2'" 


v;=i 


Both expressions indicate that 


[dX]=2- HT [dSHdL]. 




Since det(S) = det(rT'^) = Iljii it follows that 

/1 \ - -1 
[dX] = (^-j (det(S))^[dS][dL]. 

We are done. 


□ 
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Proposition 2.27. With the same notations as in Proposition 12.261 it holds that 



(2.36) 


and for m = n 




Define the normalized orthogonal measures as 


\ 2” TT 2 / \ 2" TT 2 


i>i 


(2.37) 


(2.38) 


or 

dp{V):= [dG], (2.39) 

V2"7r2 / 

where dG = V^dVfor V = [vi,...,Vn\ G 0{n). It holds that [dV"] = A(y^dy) is invariant under 
simultaneous translations V —)■ UVW, where U,W E 0{n). That is, d}i{V) is an invariant measure 
under both left and right translations, i.e. Haar measure over 0{n). 


Proof. We know that 

f [dX]e-^<^^") = 71^. 

■Jx 

From Proposition 12.261 via the transformation X = -\/SL^, we see that 

/ ldX]c- = (ty [ det(S)“F‘£-T'W|dS| X f IdL], 

Jx \2 / Js >0 JDL=tm 

We also see from the definition of Tp{ix.) that 

r,(^) = [ det(S)f-"^^e-Tr(s)[^s]_ 

V2/ Js >0 

Then 



For m = n, the result follows easily. 
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For fixed U,W G 0{n), we have 

{UVWyd{UVW) = (W^V^W) {U-dV-W) = W^ -dG- W, 

implying that 

A [{UVWyd{UVW)] = [dG]. 

That is, dpi{V) = d}i{UVW) for all U,V,W G 0(n), dpi{V) is an invariant measure under both 
left and right translations over U{n). □ 


3 Volumes of unitary groups 

3.1 Preliminary 

In Section |2l we dealt with matrices where the entries are either real constants or real variables. 
Here we consider the matrices whose entries are complex quantities. When the matrices are real, 
we will use the same notations as in Section |2l In the complex case, the matrix variable X will 
be denoted by X to indicate that the entries in X are complex variables so that the entries of 
theorems in Section will not be confused with those in Section |2l The complex conjugate of a 
matrix A will be denoted by A and the conjugate transpose by A*. The determinant of A will 
be denoted by det(A). The absolute value of a scalar a will also be denoted by \a\. The wedge 
product of differentials in X will be denoted by [dX] and the matrix of differentials by dX. 

It is assumed that the reader is familiar with the basic properties of real and complex matrices. 
Some properties of complex matrices will be listed here for convenience. 

A matrix X with complex elements can always be written as X = Xi + \/—IX 2 where 
Xi = Re(X) and X 2 = Im(X) are real matrices. Let us examine the wedge product of the 
differentials in X. In general, there are n^ real variables in Xi and another n^ real variables and 
the wedge product of fhe differentials will be denoted by the following: 


Notation 1: 


[dX] [dXi][dX 2 ] or [dX] 


d (Re 


d (Im 


(3.1) 


where [dXi] is the wedge product in dXi and [dX 2 ] is the wedge product in dX 2 . In this notation 
an empty product is interpreted as unity. That is, when the matrix X is real then X 2 is null and 
[dX] := [dXi]. If X is a hermitian matrix, then Xi is symmetric and X 2 is skew symmetric, and in 
this case 


[dXi] = and [dX 2 ] = Ay»tdxj^^ 


(3.2) 
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where Xi — 


^jk 


and X 2 = 


'{2) 

^jk 


. If y is a scalar function of X = Xi + \/—IX 2 then Y can be 


written as Y = Yi + where Yi and Y 2 are real. Thus if Y = F(X) it is a transformation 

of (Xi,X 2 ) to (Yi, Yz) or where (Yi, Yz) is written as a function of (Xi,X 2 ) then we will use the 
following notation for the Jacobian in the complex case. 


Notation 2: (Jacobians in the complex case). /(Yi^Yz : Xi,X 2 ): Jacobian of the transforma¬ 
tion where Yi and Y 2 are written as functions of Xi and X 2 or where Y = Yi -|- \/—lY 2 is a 
function of X = Xi -|- V—IX 2 . 


Lemma 3.1. Consider a matrix A G and (2n) x (2n) matrices B and C where 


A = Ai \/—IA 2 , B = 

1 

cd 

1_ 

,C = 

1- 

CM 

1 

1__ 


(N 

1 


74.2 -^1 


where Ai,A 2 G Then for det(Ai) f 0 


det(A) = |det(B)|2 = |det(C)p . 


(3.3) 


Proof. Let det(A) = a + —lb where a and b are real scalars. Then the absolute value is available 
as \l{a + \/^b){a — \/^b). If detj^i -|- \/^A 2 ) = a + y^^b, then det(Ai — '/^A 2 ) = a — 
\/—lb. Hence 


(fl -h V^b){a — y/^b) = det(Ai -|- ^/^A 2 ) det(Ai — ^/^A 2 ) 

Ai + y^A2 0 1 \ 

0 Ai — \f^^A2 J 

Adding the last n columns to the first n columns and then adding the last n rows to the first n 
rows we have 



det 


Ai -|- TA 2 0 

0 A-[ — —IA 2 


Using similar steps we have 
det 


2Ai A\ — \/-—TA 2 
Ai — V—IA2 Ai — V— IA2 


= det 


= det 


= det 


= det 


2Ai Ai — •sf — IA 2 

Ai — ^/— 1 A 2 Ai — -s/—IA2 


2Ai Ai — \/—IA 2 

—-\/^A2 \Ai — ^\/^A2 
2Ai -^^A2 
— \/^A2 |Ai 
det(Ai) det(Ai -|- A2A^^A2) 

1) = det f [ 

—A 2 Ai / \ A 2 Ai 
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by evaluating as the determinant of partitioned matrices. Thus the absolute value of det(A) is 
given by 


det(A) = det(Ai) det(Ai + A 2 A^ = |det(B)p = |det(C)|^ . 


This establishes the result. 


□ 


Remark 3.2. Now we denote Ai = Re(A) and A 2 = Im(A). Clearly both Re(y4) and Im(A) are 
real matrices. Each complex matrix A Re(A) + \/—llm(A) can be represented faithfully as a 
block-matrix 


Thus 


A 


Re(A) -Im(A) 
Im(A) Re(A) 


A* 


Re(A)’^ Im(A)'^ 

-lm(A)'' Re(A)'' 


(3.4) 


(3.5) 


Then Y = AX can be rewritten as, via block-matrix technique. 


Re(y) -Im(y) 
Im(y) Re(y) 


Re(A) -Im(A) 
Im(A) Re(A) 


Re(X) -Im(X) 
Im(X) Re(X) 


(3.6) 


From the above, we see that the mentioned representation is an injective ring homomorphism 
which is continuous, sometimes we use the following representation: 


’ Re(y) ' 


lm(y) 



Re(A) 

Im(A) 


-Im(A) 

Re(A) 


Lemma 13.11 can be reexpressed as 


Re(X) 

Im(X) 


(3.7) 


/r Re(A) -Im(A) 
\ Im(A) Re(A) 


2 


det(A) 


det(AA*) 


(3.8) 


Proposition 3.3. Let X,Y G C” be of n independent complex variables each, A G C”^” a nonsingidar 
matrix of constants. IfY = AX, then 


IfY* = X*A*, then 


[dy] = det(AA*) [dX]. 


(3.9) 


[dy*] = (-1)” det(AA*) [dXj. 


(3.10) 
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Proof. Let X = Xi + V—IX 2 , where X^ G K”, m = 1,2. Let Y = Yi + \/—VPi, where Y^ G 
]R",ra = 1,2 are real. Y = AX implies that Ym = AXm,m = 1,2 if A = A is real. This 
transformation is such that the 2 n real variables in (Yi, Y 2 ) are written as functions of the 2 m real 
variables in (Xi,X 2 ). Let 


XI = [xii,..., x„i\, XI = [Xi2, ..., X„2], 
YJ = [yii,... ,y„i], Yf = [yn, ... ,y„ 2 ]. 


Then the Jacobian is the determinant of the following matrix of partial derivatives: 

9(Yi,Yz) ^ 9(yii,...,y„i,yi2,...,i/n2) 

a(Xi, X2) d{xu, ■■■, x„i, X12, ..., X„2) ' 


Note that 


Note that 


Thus the Jacobian is 


d{Yi,Y 2 ) _ 9(yii,...,yni,yi2,--wynz) 


9(Xi,X2) 

aYi 

aYi 

aY2 

/ = det 


9 (^ 11 / • • • / ^nlf ^ 12 / • • • / ^nl) 
a(yii,...,y„i) 


= 0 = 


= A. 


3(xii,...,x„i) 

aY2 


= A, 


axi' 


A 0 
0 A 


= det(A)2. 


If A is complex, then let A = Ai + \/—lA 2 where Ai and A 2 are real. Then 

Y = Yi + ^Y2 = (Ai + V^A2)(Xi + ^X2) 
= (AiXi -A2X2) + ^(AiX2 + A2Xi) 


implies that 


Yi = AiXi - A2X2, Y2 = A1X2 +A2X1. 


Then 


aYi 

axi 

3Y2 
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Thus the Jacobian is 


( 

A.1 — A.2 

\ 


2 


7 = det ( 

A.2 -Ai 

)- 

det(A) 


det(AA*) 


which establishes the result. The second result follows by noting that [dY*] = [dYi] ( —1)” [dY 2 ] = 
(-l)'’[dY]. □ 

Proposition 3.4. Let X,Y E of mn independent complex variables each. Let A E and 

B E nonsingular matrices of constants. IfY = AXB, then 


[dY] = det(AA*) det(BB*) [dXj. 


(3.11) 


Proof Let Y = Yj + V^Y 2 and X = Xi + 7 ^X 2 . Indeed, let Y = [Yi,...,Y„] and X = 
[Xi,..., Xn], then Yj = AXj,] = 1,... ,n when Y = AX. Thus [dYJ] = det(AA*) [dXy] for each 
j, therefore, ignoring the signs, 

|dy] = nidy,] = |det(AA*)r [dXj. 

Denoting A = Ai + \/—IA 2 , the determinant is 

/ = det 

Hence the Jacobian in this case, denoting B = Bi + \/—IB 2 , is given by 

/ = det 


Ai —A2 
A2 Ai 


Bi B2 
— B2 Bi 


For establishing our result, write Y = AZ where Z = XB. That is, 

[dY] = det(AA*) ” [dZ] = det(AA*) ” det(BB*) [dX]. 
This completes the proof. 


□ 


Remark 3.5. Another approach to the fact that [dY] = det(AA*) [dZ], where Y = AX, is 
described as follows: 

^Re(Y) = Re(A)Re(Z) - lm(A)lm(Z) 


Im(Y) = Im(A)Re(Z) + Re(A)lm(Z) 


leading to 


a (^Re(Y),lm(Y)) 
a fRe(Z),lm(Z) 


Re(A)('’) -lm(A)(") 

lm(A)(") Re(A)('’) 
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where 


Re(A)W := 

' Re(A) 


, Im(A)(”) := 

' Im(A) 




Re(A) _ 



Im(A) 


Then the Jacobian of this transformation can be computed as 


/(Re(Y),Im(Y) :Re(Z),Im(Z)) 


Re(A)(") -Im(A)(”)1\ 

Im(A)(”) Re(A)(”) ) 

det(^Re(I)(”^ + ^^Im(A)(")) ^ 



That is. 


/ (^Re(Y),Im(Y) : Re(Z),Im(Z)) = det ^ = det(A) = det(AA*) 


Proposition 3.6. Let X,A,B G be lower triangular matrices where X is matrix of 
dent complex variables, A, B are nonsingular matrices of constants. Then 


Y = X + X^ 


Y = AX 


[dY] = 22''[dX], 

[dY] = 2”[dX] if the Xjfs are real] 

[d?] = 

[dY] = I ^ J [dX] if the djj's and Xjj's are real; 


Y = XB 


[d?] = (fl 

\;=i 


hi 


[d^] = n 

V;=i 


2{n-j+l) 


2{n-j)+l 


[dX], 


[dX] if the bjj's and Xjj's are real] 
Proof Results (13.141) and (I3.15I) are trivial. Indeed, note that 


Vjk = 


Ixjj, iij = k 


Xjk, it] >k 
By the definition, ignoring the sign, we have 


[d’^^j = /\j^kdyjk = A-Lidyy; /\j>kdyjk, 


(3.12) 

(3.13) 

indepen- 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 
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where di/y^ := for yjk = y^^^ + v^yj?- So for ; = 1,... ,n, we get yjj"^ = lx^^\ 

1,2. Hence the result. If xjj's are real, the result follows easily by definition. Let 

y = Yi + V^Yi,X = Xi + ^^^X^,A = Ai + V^A2,B = Bi + 

= 1 , 2 . 




m = 


where Y„„ X^, A„j, B„„ m = 1,2 are all real. 


h'/c 


Yk 


When Y = AX we have Yi = A^Xi — A 2 X 2 and Y 2 = A 1 X 2 + A 2 X 1 . The matrix of partial 
derivative of Yi with respect to Xi, that is can be seen to be a lower triangular matrix with 
flrepeated j times,; = 1,..., n, on the diagonal. Let this matrix be denoted by Gi. That is. 


aYi _ 3X2 


:=Gi = 


Ai 




Ai[i---n-l|i---M-l] 
Let G 2 be a matrix of the same structure with 's on the diagonal. Similarly, 

A 2 


axi _ 3X2 

33^ “ ^ 


:= G2 = 


A2[i|i] 


A 2 [i---n-l|i---n-l] 


Then the Jacobian matrix is given by 

3(Yi,Y2) 

3(Xi,X2) 

Let G = Gi + V—IG 2 . Then 


Gi — G 2 
G 2 Gi 


G = 


A 


A[i|l] 


A[i---n-l|i---M-l] 


where 


A = 


Ai + v^A2, a [i I i] = Ai [i I i] + V-1 A 2 [i I i],..., 

A[i • • • n^|i • • • n^] = Ai[i • • ■ n^|i • • • n^] + v^A2[i • • • n^|l ■ ■ • ?r^]. 
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Thus 


det 


Va(Xi,X2)y 


= det 


Gi -G2l\ 

G2 Gi J 


From Lemma l3Tl the determinant is available as det(G) where G = Gi + V—IG 2 . Since G is 
triangular the absolute value of the determinant is given by 


det(G) 


det(A) det(A[i|i]) 


= nd”/,! )'■ 

/=i 


det(A[i • • ■ n — l|i ■ • • n — 1]) 


This establishes (I3.16II . Another approach is presented also: Let Y = [yi,...,y„], where Yj,j = 
1, ... ,n, is the j-th column of the matrix y. Similarly for X = [Xi,...,X„]. Now Y = AX implies 
that 


Yj = AXj, j = 


That is. 


yii 


^11 


0 


0 


0 


0 

y2i 

= A 

^21 

/ 

y22 

= A 

X 22 

, . . . , 

0 

= A 

0 



. . 


yn2 


. ^”2 _ 


ynn 


^nn 


Since Y, X, A are all lower triangular, it follows that 


yii 

y2i 

= A 

^11 

^21 

/ 

yii 

= A[i|i] 

X 22 

y^i 


^nl 


. Vnl _ 


. _ 


ynn = A[l...n-l\l...n-l]x„n, 


where A[ij\iJ\ stands for a matrix obtained from deleting the i,/-th rows and columns of A, 
respectively We can now draw the conclusion that 



det(A[i...;-l]|l...;-l])A[l...;-l]|i...;-l])* 


[dXy], 
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that indicates that 


[dY] = n[d^;]=n det(A[l.../-l]|i...;-l])A[i...;-l]|l...;-l])* [dX^-] 


det(AA*) det(A[l|l]A[i|i]^ 
2 


\an„al„\ [dX] 


— |flllfl22 ■ • • Ann I X | fl22fl33 ' ' ’ X ■ ■ ■ X [dX] 


= (ni?//P'jidxi- 


If the Xjj's and ajj's are real then note that the Xjk's for j > k contribute ajj twice that is, corre¬ 
sponding to and x^^\ whereas the ajj's appear only once corresponding to the since the 
are zeros. This establishes (I3.17D . If Y = XB and if a matrix Hi is defined corresponding to 
Gi then note that the appear n — ; + 1 times on the diagonal for j = 1,... ,n. Results (I3.18D 
and (I3.19D are established by using similar steps as in the case of (I3.16I) and (I3.17D . □ 


Proposition 3.7. Let X G be hermitian matrix of independent complex entries and A G be a 
nonsingular matrix of constants. IfY = AX A*, then 


[dY] = det(AA*) 


[dX], 


(3.20) 


Proof. Since A is nonsingular it can be written as a product of elementary matrices. Let Ei,..., Ej. 
be elementary matrices such that 

A = E,E,_i... El A* = ... E*,. 

For example, let £i be such that the j-th row of an identity matrix is multiplied by a scalar 
E = a + ^/—Ib where a,b G IR. Then EiXE^ means that the j-th row of X is multiplied by 
a -h \/—lb and the j-th column of X is multiplied by fl — \/—lb. Let 

Hi = EiXEi*, H2 = £2 HiE|, ..., = Efca^_lE^ 

Then the Jacobian of Y written as a function X is given by 

/(Y:X) = /(Y:H,_i).../(Hi:X). 

Let us evaluate [dHi] in terms of [dX] by direct computation. Since X is hermitian its diagonal 
elements are real and the elements above the leading diagonal are the complex conjugates of 
those below the leading diagonal, and Hi is also of the same structure as X. Let Ui = U-\- \/—lV 
and X = Z + where U = [ujk], V = [vjk], Z = [zjk], W = [iVjk] are all real and the diagonal 

elements of V and W are zeros. Take the ujj's and Zjj's separately. The matrix of partial derivatives 
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of wii,..., u„„ with respect to Zn,..., z„„ is a diagonal matrix with the ;-th element and 

all other elements unities. That is. 


1 


3(diag(U)) _ 9 (mii,..., ^nn) 
a(diag(Z)) 9(zii,... r Z)jn) 


_|_ ^2 _ |^|2 


:= C, 


1 


where diag(X) means the diagonal matrix, obtained by keeping the diagonal entries of X and 
ignoring the off-diagonal entries. 

The remaining variables produce a ”^”2 ^ matrix of the following type 


d(Uo, Vq) _ Aq Bq 
9(Zo,Wo) — Bq Aq 


where Uq, Vq, Zq, Wq mean that the diagonal elements are deleted, Aq is a diagonal matrix with 
n — 1 of the diagonal elements equal to a and the remaining unities and Bq is a diagonal matrix 
such that corresponding to every a in Aq there is a fc or — fc with / — 1 of them equal to —b and 
n — i oi them equal to b. Thus the Jacobian of this transformation is: 


/(Lfi : X) 


9(diag(U),Uo,Vo) 

9(diag(Z),Zo,Wo) 



/ 

’ c 

0 

0 

\ 

det 


0 

^0 

So 



V 

0 

-So 

Aq 

/ 


det(C) det 


^0 So \ 
— Bo ^0 _ / 


From Lemma 13.11 the determinant is 


det(Ao -|- ^/^Bo) 1^. That is. 


det 


/9(LJo,Vo)\ 

V9(Zo,Wo)y' 


det 


^0 So \ 

— So Aq _ ) 


= det((Ao + ^So)(Ao + y^So)*) 
= {a^ + b^)"-\ 


Thus 


[dlii] = {a^ + b^)^[dX] 


n 


det(£i£j) 


[dX]. 


Note that interchanges of rows and columns can produce only a change in the sign in the deter¬ 
minant, the addition of a row (column) to another row (column) does not change the determinant 
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n 


and elementary matrices of the type Ei will produce 
puting /(til : X),J{U 2 ■ Ui) etc we have 


det(EiE*) 


[dY] 


det(AA*) 


n 


[dX], 


in the Jacobian. Thus by com- 


As a specific example, the configuration of the partial derivatives for n = 3 with J = 2 is the 
following: If 


then 



Zll 

* 

* 


1 0 

0 

X = 

Z21 + 3 /^ 71)21 Z 22 

* 

and El = 

0 c 

0 


Z3I + Z32 + 3/^W32 

2:33 _ 


0 0 

1 


Zll 


* 

* 


Ui = 

C(Z21 + V^W23) 


|c|^Z22 

* 



Z31 + V — 1 h ^31 

C(Z 32 + \/^W32) 

Z33 



l/ll 

* 

* 




M 2 I + a/-1^21 M 22 * 

W3I + M32 + \/^V 32 W33 


thus 


Now 


Now 


Wll — 2:11,1/22 — |c|^Z 22 , W33 — 2:33, 
U 21 = aZ 21 - bW21,U31 = Z31,W32 = aZ32 + bW32, 


V 21 = 

azv 2 i + bz 

21/ 

^^31 

— 

“^31/J^32 


aW32 — 

bZ32. 

Mil 


1 

0 

0 

0 

0 

0 

0 

0 

0 ' 


Zll 

U 22 


0 

\cf' 

0 

0 

0 

0 

0 

0 

0 


Z 22 

“33 


0 

0 

1 

0 

0 

0 

0 

0 

0 


Z 33 

“21 


0 

0 

0 

a 

0 

0 

-b 

0 

0 


Z 21 

“31 

= 

0 

0 

0 

0 

1 

0 

0 

0 

0 


Z 31 

“32 


0 

0 

0 

0 

0 

a 

0 

0 

b 


Z 32 

U 21 


0 

0 

0 

b 

0 

0 

a 

0 

0 


W 21 

i^31 


0 

0 

0 

0 

0 

0 

0 

1 

0 


W 31 

V 32 


0 

0 

0 

0 

0 

—b 

0 

0 

a 


W 32 


C = 

1 

|c|^ 

/ Ao = 

a 

1 

, Eq = 

1 

1 

0 

_ 1 


1 


a 


b 
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We are done. 


□ 


Remark 3.8. If X is skew hermitian then the diagonal elements are purely imaginary, that is, the 
real parts are zeros. It is easy to note that the structure of the Jacobian matrix for a transformation 
of the type Y = AXA*, where X* = —X, remains the same as that in the hermitian case of 
Proposition 13.71 The roles of {ujj,Zjj)'s and are interchanged. Thus the next theorem 

will be stated without proof. 


Proposition 3.9. Let X G skew hermitian matrix of independent complex entries. Let A G 
be a nonsingular matrix of constants. IfY = AXA*, then 


[dY] = det(AA*) [dX]. 


(3.21) 


Some simple nonlinear transformations will be considered here. These are transformations 
which become linear transformations in the differentials so that the Jacobian of the original trans¬ 
formation becomes the Jacobian of the linear transformation where the matrices of differentials 
are treated as the new variables and everything else as constants. 

Proposition 3.10. Let X G fee hermitian positive definite matrix of independent complex variables. 
Let T G be lower triangular and Q G be upper triangular matrices of independent complex 
variables with real and positive diagonal elements. Then 

X = ff* [dX] = 2” [df], (3.22) 

X = QQ* [dX] = 2” [dQ]. (3.23) 


Proof. When the diagonal elements of the triangular matrices are real and positive there exist 
unique representations X = TT* and X = QQ*. Let X = Xi -|- V—IX 2 and T = Ti + V—IT 2 , 
where X = [xjj,],f = \ijk],fk = 0,j < k,Xm = = 1,2. Note that Xi is 

symmetric and X 2 is skew symmetric. The diagonal elements of X 2 and T 2 are zeros. Hence when 
considering the Jacobian we should take Xjj, j = 1,... ,p and xyj., j > k separately. 

X = ff* Xi + ^/^X2 = (Ti + V^T2){TJ - ^^TJ) 

^ ( Xi = TiTJ + T 2 TI 

\X2 = T2T{ - TiTl 


with t^j'^ = tjj, =0,/ = l,...,n. Note that 


= f (t]l^ 


Aitf\ 


+ 1 ( i,i 


( 2 ) 


+ ( 
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implies that 


So 






( 1 ) 


= =2tyy,/ = l,.. 


n. 




2fii 


:= Z. 


2t 


nn 


Now consider the 


s for j > k. It is easy to note that 


3(Xio,X2o) _ U V 
d{Tio,T2o) W Y 


where a zero indicates that the x^j'^'s are removed and the derivatives are taken with respect to the 
fj^^'s and for j > k. U and Y are lower triangular matrices with tjj repeated n — j times along 
the diagonal and V is of the same form as U but with = 0 along the diagonal and the 
replaced by the For example, take the xj^^'s in the order x^^,x^j^,..., ,x^^,.. 

and also in the same order. Then we get the ”^”2 ^ matrix 

fii 0 ■ ■ ■ 0 

* fii • • • 0 

: : 0 

* * ■ • ■ tn-l,n-l 

where the *'s indicate the presence of elements some of which may be zeros. Since U and V are 
lower triangular with the diagonal elements of V being zeros, one can make W null by adding 
suitable combinations of the rows of {U,V). This will not alter the lower triangular nature or the 
diagonal elements of Y. Then the determinant is given by 


U = 


3Xio 
3Ti 


10 


det 

W Y 

Multiply with the 2tjj's for ; = 1,..., n to establish the result. As a specific example, we consider 
the case where n = 3. Let X G Denote X = Xi + ^/^X 2 . Thus 



1 

y(l) 

X21 

-^31 


1 

0 

y(2) 

X21 

r(2) 1 

■^31 

II 

^( 1 ) 

X21 

r(l) 

X22 

y(l) 

-^32 

, X2 = 

J2) 

X21 

0 

r(2) 

■*•32 


JY 

■^31 

^( 1 ) 

■^32 

^( 1 ) 

■^33 


1 

J2) 

Hi 

0 


= det(D) det(Y) =n^/) 


2(n-;) 


/=1 
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Similarly, let T = Ti + \/—lT 2 . We also have: 



hi 

0 

0 


0 

0 

0 

Ti = 

hi 

h2 

0 

, T2 = 

hi 

0 

0 


L 31 

h2 

h3 J 


L 31 

h2 

0 


Now Xi = TiT^ + T 2 TI can be expanded as follows: 


y(i) 

^(1) _ 
X 22 — 

r(l) _ 

-^33 “ 


- f 


.(1) 

11 


dl) _ dl)dl) 


'■21 


— ‘■21 hi ' -^31 


(1) _ dl)dl). 

- /-Ol ti 


*■31 ' 


“ y-ii j ' \^h 2 j ^ l^hi j ' -^32 ~ hi hi ^ h 2 ‘■22 ' hi hi' 

= (410 




.(1) 

'33 


+ t 


.( 2 ) 

'31 


+ t 


.( 2 ) 

'32 


Then X 2 = T 2 TJ — T 1 T 2 can be expanded as follows: 

y(2) _ h2),(l) (2) _ (2) (1) (2) _ (2) (1) (2) (1) 
•^21 “ hi hi ' -^31 “ 31 11 ' -^32 ~ hi 21 ' h2 h2 ' 

From the above, we see that 


dY^j^ 


/•hi 

0 

0 

0 

0 

0 

0 

0 

0 


dY^^ 


0 

Z.I 22 

0 

Zhi 

0 

0 

Zhi 

0 

0 




0 

0 

^h3 

0 

^hi 

2 f*'^^ 

Zh2 

0 

2 f^^^ 

Zhl 

2 f*'^^ 

AI 32 


dY2^^ 


hi 

0 

0 

hi 

0 

0 

0 

0 

0 


dY^P 

= 

hi 

0 

0 

0 

hi 

0 

0 

0 

0 


dY^^ 


0 

h2 

0 

hi 

hi 

h2 

hi 

hi 

0 


dY2^^ 


hi 

0 

0 

0 

0 

0 

hi 

0 

0 


d4? 


hi 

0 

0 

0 

0 

0 

0 

hi 

0 


. ^^32 _ 


0 

h2 

0 

hi 

0 

0 

0 

hi 

h2 J 



hi 

d41' 

h3 

hi 

Uhl 

h2 
hi 

Uhl 


dt 


d4 

d4 


d4 

d4 


df 


( 2 ) 

32 
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In what follows, we compute the Jacobian of fhis transformation: 


J(X : T) 


/ 


0 

0 

0 

0 

0 

0 

0 

0 

\ 


0 

^‘22 

0 

^‘21 

0 

0 

^‘21 

0 

0 



0 

0 

^‘33 

0 

^‘31 

^‘32 

0 

^‘31 

Z‘32 



‘21 

0 

0 

‘11 

0 

0 

0 

0 

0 



‘31 

0 

0 

0 

‘ll 

0 

0 

0 

0 



0 

‘32 

0 

‘31 

‘21 

‘22 

‘31 

‘21 

0 



‘21 

0 

0 

0 

0 

0 

‘11 

0 

0 



‘31 

0 

0 

0 

0 

0 

0 

‘ll 

0 


V 

0 

‘32 

0 

‘31 

0 

0 

0 

‘21 

‘22 J 

/ 


/ 

1 

0 

0 

0 

0 

0 

0 

0 

0 

\ 


0 

2f^^^ 

^‘22 

0 

^‘21 

0 

0 

2f^^^ 

^‘21 

0 

0 



0 

0 

2f^^^ 

^‘33 

0 

^‘31 

^‘32 

0 

2f^^^ 

^‘31 

2f^^^ 

^‘32 



‘21 

0 

0 

‘ll 

0 

0 

0 

0 

0 



‘31 

0 

0 

0 

‘ll 

0 

0 

0 

0 



0 

‘32 

0 

‘31 

‘21 

‘22 

‘31 

‘21 

0 



‘21 

0 

0 

0 

0 

0 

‘ll 

0 

0 



‘31 

0 

0 

0 

0 

0 

0 

‘ll 

0 


V 

0 

‘32 

0 

‘31 

0 

0 

0 

‘21 

‘22 J 

/ 


by adding the corresponding multiples of the first row to the second row through the last one, 
respectively, we get 



/ 

’ 1 

0 

0 

0 

0 

0 

0 

0 

0 

\ 



0 

Z‘22 

0 

2f^^^ 

^‘21 

0 

0 

^‘21 

0 

0 




0 

0 

^‘33 

0 

^‘31 

^‘32 

0 

^‘31 

^‘32 




0 

0 

0 

‘ll 

0 

0 

0 

0 

0 


21^^ det 


0 

0 

0 

0 

‘ll 

0 

0 

0 

0 




0 

‘32 

0 

‘31 

‘21 

‘22 

‘31 

‘21 

0 




0 

0 

0 

0 

0 

0 

‘11 

0 

0 




0 

0 

0 

0 

0 

0 

0 

‘ll 

0 



1 

0 

‘32 

0 

‘31 

0 

0 

0 

‘21 

‘22 J 

/ 


iteratively, finally we get the final result. We also take a simple approach (i.e. by definition) with 
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a tedious computation as follows: 


and 


dx22^ 

dxli^ 

.4' 


2fiidfii, dx2^^ = f2i^dtii + fiid4i^ dxgp = tg^^dfn + fndlg^; 

2f22df22 + 2 t^2idt^2i + 245^d4?, 

44^^44 + 44^44 + 44^^22 + f22d44 + 44^44 + 44 ^ 44 ' 
244 d44 + 244 d44 + 2f33df33 + 244 d44 + 244 d44' 


dx^^ = 44dlii + liid44' dxg^^ = 44dfii + fiid44/ 
d^34 ~ 44d44 + 44d44 + ^22d44 + 44di22- 


Hence we can compute the Jacobian by definition as follows: 


[dX] = dx^j^ A dx)^-^ A dXg^'' A dx^^-* A dXg^'' A dx^^"* A dxg^'' A dxg^'' A dxgg 


.(1) 


.(1) 


.(2) 


.(2) 


.(1) 


.(2) 


.(1) 


( 1 ) 


23 n [df|. 


V /-1 


The proof in the case of X = QQ* is similar but in this case it can be seen that the triangular 
matrices corresponding to U and Y will have tjj repeated 7 — 1 times along the diagonal for 
7 = l,...,n. □ 


Example 3.11. If X* = X G is positive definite , and Re(a) > n — 1, 


r„(a) 


[dX] 


det X 


-Tr(x) 


JX>0 ^ / 

7r“^^r(a)r(a — l)---r(a — n + 1 ). 


(3.24) 


Indeed, let T = [Fyjt]/ tjk = 0,j < k he a lower triangular matrix with real and positive diagonal 
elements tjj > 0,j = 1,... ,n such that X = TT*. Then from Proposition |3T0] 

[dX] =2” [df] 




Note that 


Tr (x) =Tr (rr) = + \t 2 i\'+ ■ ■ ■ + ■ ■ ■ + 


M 


and 


det(x)” [dx ]=r (n 4 r 


d=i 
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The integral over X splits into n integrals over the tjj's and ”^”2 integrals over the j > k. 
Note that 0 < tjj < 00 , — 00 < < 00 , — 00 < < 00 , where But 


^ 2 j+i^ = Y{a-i + 1), Re(a) > j-I, 


for j = 1,... ,n, so Re(a) > n — 1 and 


/* I—- |2 pOO pOO 






The desired result is obtained. 


Definition 3.12 (T„(a): complex matrix-variate gamma). It is defined as stated in Example 13.111 
We will write with a tilde over T to distinguish it from the matrix-variate gamma in the real case. 


Example 3.13. Show that 


/W = 


det(B) det(X) 


.Tr(BX) 


T„{a) 


for B* = B > 0,X* = X>0, Re(a) > n — 1 and /(X) = 0 elsewhere, is a density function for X 
where B is a constant matrix, with r„(a) as given in Definition 13. 121 Indeed, evidently /(X) ^ 0 
for all X and for all X and it remains to show that the total integral is unity. Since B is hermitian 
positive definite there exists a nonsmgular C such that B = C*C. Then 


Tr BX =Tr CXC 


Hence from Proposition l3.10l 

y = CXC* [dY] = 


det CC 


[dX] = det(B) [dX], 


and 


X = C^WC* 


det(X) =|det(CC*)r^ det(y) 


Then 


ix>o 


r r 

det(y) 

i [dY]- 
/y>o 

T 


-Tr(Y) 


Tn{a) 


= 1 . 


But from Example 13.1 li the right side is unity for Re (a) >n — 1. This density /(X) is known as 
the complex matrix-variate density with the parameters a and B. 
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3.2 The computation of volumes 


Definition 3.14 (Semiuniatry and unitary matrices). A p x n matrix U is said to be semiunitary 
if UU* = Ip for p < n or U*U = 1„ for p > n. When n = p and UU* = then U is called a 
unitary matrix. The set of all n x n unitary matrices is denoted by W(n). That is. 


U{n) := ju G : UU* = !„}, 


(3.25) 


where denotes the set of all n x n complex matrices. 


Definition 3.15 (A hermitian or a skew hermitian matrix). Let A G If A = A*, then A is 

said to be hermitian and if A* = — A, then it is skew hermitian. 


When dealing with unitary matrices a basic property to be noted is the following: 


UU* = 1 ^ U*dU = -dU*U. 


But \^U*dUj = dU*U, which means that U*dU is a skew hermitian matrix. The wedge product 

of U*dU, namely, A (jU*dU^ enters into the picture when evaluating the Jacobians involving 
unitary transformations. Hence this will be denoted by dG for convenience. Starting from 
U*U = 1„ onehas dG-H*. 


Assume that U = [uy] G U (n) where uy G C. Let Ujj = | Ujj | by Euler's formula, where 

9j G [—n,n]. Then 




U12S • 

■ ui„e 


■ gv/^01 

0 

0 

u = 


1M 22 1 



0 

(>\AU02 

0 


Unie^'^^'^ 


1 W)2n 1 


0 

0 



(3.26) 


This indicates that any U G U{n) can be factorized into a product of a unitary matrix with 
diagonal entries being nonnegative and a diagonal unitary matrix. It is easily seen that such 
factorization of a given unitary matrix is unique. In fact, we have a correspondence which is 
one-to-one: 


U{n) ~ {U(n)IU{\Y^) X U{\Y^. 


Notation. When LI is a n x n unitary matrix of independent complex entries, U* its conjugate 
transpose and dii the matrix of differentials then the wedge product in dG := dlJ ■ U* will be 


denoted by [dG]. That is, ignoring the 

sign. 




< 

il 

'dU-U*"^ 

= A 

'u-dU*" 

. 


44 









If the diagonal entries or the entries in one row of this unitary matrix U are assumed to be real, 
then the skew hermitian matrix dlJ ■ U* will be denoted by dGi and its wedge product by 


[dGi] := A (did - G*) . 


Indeed, [dG] here means the wedge product over U{n), but however [dGi] means the wedge 
product over Zd(n)/W(l)^’k Therefore, for any measurable function/over Zd(n), [dG] = [dGi][dD], 


f f{U)[dG]= [ [ 

Ju{n) JUi{n) Ju{l)x’' 


/(yD)[dGi][dD], 


where U = VD for V G Ui{n), dG = U*dU and dGi = V*dV. Furthermore, let V = [vij] for 
Vij G C and Vjj = Vjj G R+. Then it holds still that V G U{n), thus Vjj is not an independent 


variable, for example, Un = y/1 — j ^21 j —••• — [ Vni \ ■ From this, we see that 


[dGi] = (Re(y*dy)y) d (lm(y*dy); 

i<i 


and 

[dG] = ( fjlm(;j*dJi)y^ I X ]QRe(Ji*d!J)yIm(Ji*d!J)y. 

/ i<i 

Proposition 3.16. Let T G fee loiver triangular and U ^U{n) be of independent complex variables. 
Let X = fU. Then: 


(i) for all the diagonal entries tjj,] = 1,..., n 0 / T being real and positive, 

[dx] = [df][dG] 

where dG = dU ■ U*; and 

(ii) for all the diagonal entries in U being real, 

[dX]= ■[df][dGi] 

zvhere dGi = dU ■ U*. 

Proof. Taking differentials in X = TU one has 

dX = df ■ Lf + f ■ dQ. 

Postmultiplying by U* and observing that UU* = we have 

dX ■ Lf* = df + f ■ dLf ■ LI*. 


(3.27) 


(3.28) 


(3.29) 
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(i). Let the diagonal elements in T be real and positive and all other elements in T and U be 
complex. Let 


dV = dX-U* [dy] = [dX] (3.30) 

ignoring the sign, since U is unitary. Let dG = dli ■ U* and its wedge product be [dG]. Then 

dy = df + f-dG (3.31) 

where dG is skew hermitian. Write 


y = 
f = 


[Vjkh Vjk = 


h2) 


;; =tn>0. 


dG = [dgjk], dgjk = dg]l^ + \/^dg^l\ dgf^^ = -dgl^\ dg^l^ = dgl^\ 

From Eq. (|3.31ll . 

I dtj^- + (^tjidgik + • • • + tjjdgjk) , i ^ k 
■+^jidgjk), i<k. 


dvjk = 


1 (Jn^gik 


The new variables are dyj^^ = = g^’i^\i ^ k,m = 1,2. The matrices of 

partial derivatives are easily seen to be the following: 


dTjj 


= 1 , 


> k 


^k 

^^jk 


,j > k 


= 1, m = 1,2, 


= A, 


dv 


( 2 ) 




,j > k 


= B 


where A and B are triangular matrices with tjj repeated n — j times, 

= diag(fii,.. .,f„n). 




By using the above identity matrices one can wipe out other submatrices in the same rows and 
columns and using the triangular blocks one can wipe out other blocks below it when evaluating 
the determinant of the Jacobian matrix and finally the determinant in absolute value reduces to 
the form 


det(A) det(B)fi 


f — I I f2(«“;)+i 
1 ' ' ' ‘■nn — J_ ‘■jj 
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Hence the result. As a specific example, we consider the case where n = 3. We expand the 
expression: dV = dT + T ■ dG. That is. 


Furthermore, 


and 


Thus 



dW 

= 


dy2 

= 

dt;d) 

0^13 

dv^^2 

dyg^ 

dt7d) 

0^23 

dV22 

dZ7^3^ 

dyd) 

UC 33 

^^31 

dt^33^ 


— T 2 d (^Im 
+ Tid fim 


•4.4 

1_ 

0 

1 

0 


1 

0 

1 

0 


1 

0 



UC21 

dt^2^ 

0 

+ 

''21 

fd) 

'22 

0 


-dgii^ 

0 

dg23^ 

L ^'^31 

dt^2^ 

d43^ . 


fd) 

L 31 

fd) 

'•32 

fd) 
''33 J 


. -d^is^ 

-dg23^ 

0 


-1 

0 

0 

1 

0 


^Sn 

^sfi 

d^S? 

fd) 

'21 

0 

0 


d§u 

^g?2 

d^2? 

fd) 

L 31 

fd) 

'32 

0 


_ d^B 

^gf3 

d4? . 


dyd) 

UCii 

dv^ 

dz;d^ 





d^? 

dv^ 

dug) 





dyd) 

dyd) 

dz^3? 





0 

0 

0 


0 

0 

0 

dfd) 

0121 

0 

0 

+ 

fd) 

'-21 

0 

0 

dfd^ 

0131 

d4? 

0 


fd) 

'^31 

fd) 

'32 

0 


0 

0 


d^B 

d^23^ 



r fd) 

'■11 

0 

0 


d^n 

d^S? 

^S?3 

+ 

fd) 

'21 

fd) 

'22 

0 


d^!? 

d^g^ 

^§S 


fd) 

L 31 

fd) 

'32 

OJ-^ 

1_ 


. dS?3 

d^g^ 

dg3? . 


d^Sl^ = 
d^ 2 l^ = 

d^23^ = 

d4i> = 

-g- = 


dfii, dr;{ 2 ^ = fndgi^^ du^g^ = hid^Jg^ 

d41^ - f22dg{2^ - 4?dg{5^ d42 ^ df22 + 44d^S2^ - 4?d^52^ 

+ ^22dg^3^ - 4?dgj3\ 

AfW _ fd) jJl) _ X Jad) _ fd) j J2) _ .(2) 1 J2) 
dtsi r32 0^32 *^330^33 133 0^33 1^2 ^15^12 ' 

Hfd) , xd) jdl) _ X J Jl) _ x(2) J J2) _ x(2) jd2) 

0132 + ^33 0^32 1330^23 ^31 Clg32 132 0^22 ' 


rl7;d) _ ^fd)jdl) , 3(1 )jJ1 )_3(2) n f2 )_,(2) , (2) 

d^33 — dr33 + r33 0^33 + t^2 ^§23 hi ds^i3 hi ds^23 
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and 


dvf^ = hidgf^, = tudg^^\ dvf^ = hidgf^, 
dv^S = d4? + 41^dgii^ + f22dgj2^-4?dg{2^ 
dy22^ = + 44^^12^ + t22dg^^\ 

d44 - 44dg!4+44dgi?+ ^22dgg\ 

dvf^ = d44 - 4? dgl4 + 44 dg{4 + 44 dgl4 + ^ssd^iJ, 

d^? ^ d44 + 44dgl4 + 44dg{4 + 44 d^S^ + ^ssd^g^ 

d^? = 44dg{4+44 dg24 + 44dg{4 + + ^33dg54- 

According to the definition, we now have dvjk = dyj^^ A dv^f^, then 

[dV] = drill A drii2 A d^'is A dr;2i A d522 A dc '23 A drisi A d532 A dvo,^, 

= IdflldG]. 

(ii). Let the diagonal elements of U be real and all other elements in U and T complex. Starting 
from Eq. (I3.31I) , observing that dG is dGi in this case with the wedge product [dGi], and taking 
the variables dvjj^'s in the order dv^^\i ^ k,dv^^\i < Kdvf,i ^>dvf,i < k and the other 

variables in the order dt^l\j ^ k,dg^l\i > k,dt^l\j ^ k,dg^f^,j > k, we have the following 

configuration in the Jacobian matrix: 

1 * * * 

0 —Ai 0 A 2 

0 * 1 * 

0 —A 2 0 —Ai 

where the matrices marked by * can be made null by operating with the first and third column 
submatrices when taking the determinant. Thus they can be taken as null matrices, and Ai and 
A 2 are triangular matrices with respectively and repeated n — j times in the diagonal. The 
Jacobian matrix can be reduced to the form 


'a b ' 

,A = 

' 1 

0 


-B A 


0 

—Ai 



0 0 
0 A 2 


Then the determinant is given by 


A B 
-B A 


det((A + V^B){A + V-IB)*) 


det((-Ai + ^/^A2){-Al + ^/^A2) 

niiT-"-'' 

/=! 
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since — Ai + \/—XA^ is triangular with the diagonal elements + V—repeated n — j 

times, giving repeated n— j times in the final determinant and hence 

the result. □ 

By using Proposition 13.161 one can obtain expressions for the integral over U of [dG] and 
[dGi]. These will be stated as corollaries here and the proofs will be given after stating both the 
corollaries. 

Theorem 3.17. Let dG = dLf • U*, where U ^U{n). Then 



Proof. Let Xhe an x n matrix of independent complex variables. Let 

.7 X J X 


= n 


smce 


JXjl 


e iLd dxj^ = 


' —CO J — 




= TC. 


Consider the transformation used in Proposition 13.161 with tjj's real and positive. Then 

Tr (XX*] =Tr( ff*" 


and let 


B = 




But 


and for j > k 


pOO -N -1 2 1 

/ f : ” ^ ^ e~^i>dtjj = -Y{n — j + 1) for n — ; + 1 > 0 
Jo " 2 


/ 

Jtn 


e dfyjt = TT. 


Then the integral over T gives 


2-v'^nr(^-; + l) =2-'T„(n). 
M 
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Hence 


r _ 2"7r” 

/ [dG] = 

Iu(n) r„(n) 


□ 


Theorem 3.18. Let Hi G U{n) with the diagonal elements real. Let dGi = dJii ■ U^. Let the full unitary 
group of such n X n matrices Ui be denoted byUi{n) = U{n) Then 


vol 


(Hi(n)) =Yo\{U{n)/U{\Y^) = f [dGi] = ^ 


n{n—l) 


n(n-i) 

n 2 


fn) l!2 !---(n-l)!' 


(3.33) 


Proof. Now consider the transformation used in Proposition 13. 161 with all the elements in T com¬ 
plex and the diagonal elements of U real. Then 


B = 




Note that 


nfi 

j>k 


n(n-l) 


e IHI dfyjt = TT 2 


Let tjj = t = ti + ^/—lt 2 . Put fi = rcos0 and t 2 = rsind. Let the integral over tjj be denoted by 


aj. Then 


a; = 


= 4 


e-\^\ dJ = 

+ CX) /* + oo 


= 4 


/ -foo /*+oo ... 

-00 J —oo 

pd-00 /‘ + 00 

Jo Jo 

{r^Y~! ■ e-’"' • r ■ drdO 


dti 


Je=:0 Jr=0 
TTT{n — i + 1) for n — i + 1 > 0. 


Then 

implies that 

which establishes the result. 
Example 3.19. Evaluate the integral 


.(«) /- [dG 

JUi 


ij = 71" 


L [dGi 
Ju, 


lUi T„{n) 


□ 
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A(a) = 


= hdX] det{XX* 

Jx 


g-Tr(xx*) ^ + n) 

fn{n) 


(3.34) 


for Re(a) > —1, where X G matrix of independent complex variables. Indeed, put X = TU, 
where U G U{n), T is lower triangular with real distinct and positive diagonal elements. Then 


det(XX*; 
Tr (XX* 


fit' 

j=i 




;>/c 


M i>^ 


[dX] 

/jdG] 

Ju 


r^' 


Thus 


A(a) = y^_[dT][dG] 


' n po 

Bi 


j2(i+2(n-;)+lg 


n/ 


e \ id dt 


jk 


'U 


.[dG] 


^ \ «("-!) 2”7r" 

= 2 "nr(a + n-/ + l) -TT 2 

V 7=1 J r„(n) 

for Re(a + n) > n — 1 or Re(a) > —1. That is, 

A(a) = 77 -”^for Re(a) > —1. 

r„(pn) 

Proposition 3.20. Let X G fl hermitian matrix of independent complex entries with real distinct 

eigenvalues Ai > Ai > ■ ■ ■ > A„. Let U ElAin) with real diagonal entries and let X = UDU*, where 
D = diag(Ai,..., A„). Then 


[dX] = (ny>;c|A^-A;f)-[dD][dGi], 


(3.35) 


where dGi = U* ■ dU. 


Proof. Take the differentials in X = UDU* to get 

dX = dQ • D • Lt* + Lf • dD • Lt* + Lt • D ■ dLT*. 
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Premultiply by U*, postmultiply by U and observe that dGi is skew hermitian. Then one has 

dW = dGi • D + dD - D ■ dGi 

where dW = U* ■ dXU with [dW] = [dX], Using the same steps as in the proof of Proposi¬ 
tion |5T0l we have 

[dW] = [dD][dGi]. 

Hence the result follows. □ 

Example 3.21. Let D = diag(Ai,..., A„) where the Ay's are real distinct and positive or let Ai > 
■ ■ • > A„ > 0. Show that 


(i) 


/Ai>--->A„>0 


|dD](n|A.-A,rp-^'<°) = A^=(nrO)| ; 


d>k 




(ii) 


[ [dD](Yl\\k 


- A; 


Vi=i 


-n\ _ -Tr(D) _ r„(tt)r„(n) 


In fact, let y be a n X n hermitian positive definite matrix, U a unitary matrix with real diagonal 
elements such that 


U*YU = D = diag(Ai,...,A„). 

From the matrix-variate gamma integral 

— L _ [dY] det(y) for Re(ft;) > n — 1. 

Jy=y*>o 


Hence 


Jy=Y*>0 


Then Y = UDU* implies that 


[dY] = i Ale - Ay 1^ j [dD] [dGi], dGi = U* ■ dU. 
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Note that 


Tr y = Tr UDU* =Tr(D), 


det(Y) = Ini'^d 


Then 


T„{n) = [ [dD] x / [dGi], 

7Ai>->A„>0 " ) Ju 


Clearly 




TT‘ 


n{n—l) 


Tnin) 


Substituting this, results (i) and (ii) follow. 

Remark 3.22. We can try to compute the following integral in (I3.35I) : 

L [dx] = [ [ [dG^] 

JX:Tr{x)=l ./Ai>A2>-->A„>0 ^ J t<) f=l JUi{n) 

which is equivalent to the following 




1 npor(n-;)r(n-/ + i) 


n! 


r(n2 


nr.,ro) 


1 .1.-11 rnffC2Y--rfM + n 

n!^ r(n2) 


where we used the integral formula: 


;=1 / i<j ;=1 


n;to'r(n-;)r(«-; + i) 

r(n2) 


This means that 


vol (D (C”)) = n 


;^ T(1)T(2)---T(n) 

T(n2) 


(3.36) 
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We make some remarks here: to obtain the volume of the set of mixed states acting on C”, 
one has to integrate the volume element [dX], By definition, the first integral gives where 

Cjj*® = and 

n 

1 poo poo ( n \ yi _ ft 

= I ■■■I ^ nArn(A--A;)^ndA; 

Q V ;=1 / k i<j 7=1 

1 " r(i+,-f)r (» + (/-i)f) 

r («« + j=i + 

while the second is equal to the volume of the flag manifold. To make the diagonalization 
unique, one has to restrict to a certain order of eigenvalues, say Ai > A 2 > • • • > A„ (a generic 
density matrix is not degenerate), which corresponds to a choice of a certain Weyl chamber of the 
eigenvalue simplex In other words, different permutations of the vector of n generically 

different permutations (Weyl chambers) equals to n!. This is why the factor ^ appears in the 
right hand side in the above identity. In summary, the transformation 

(d,Q) 


such that X = UDU*, is one-to-one if and only if D = diag(Ai,..., A„), where Ai > • • ■ > A,, 
and U e U{n)IU{\Y^. 

Remark 3.23. In this remark, we will discuss the connection between two integrals HI: for a, /3 > 
0 , 

n 

' foo 7 m / n \ n n 

= L "L ■’h-E-'/ n-'r' n (A,-A,)PndA/, 

\ ;=1 / k=l 7=1 

n 

. . ' /-OO 7 m ( n \ n n 

X^Y\oc,^) - / expi-f^Xj\Yl4~^ n 

\ 7=1 / lc=l lsSi<7sSn 7=1 

We introduce an auxiliary variable t in the expression of /3), and define I{t) as 

n 

' nco 7^ / n \ n n 

L m^-ea, riAr* n (A,-A,)/*ndA,. 

•A'A ■A'J \ ;=1 / A:=l l<i</<n /=1 


(3.37) 
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Then 7(1) = Taking the Laplace transform, denoted by Zf, of 1(f), we obtain 

n 

poo poo p / ft \ n _ 

if(l) := / ■■■/ n (A,-A,)^ndA; (3.38) 

I"' V 1=1 / J k=l 7=1 

n 

' roo ^ ^ / n \ n n 

= L " L «p n 

\ /=1 / i:=l 1<!<7^« 7=1 


= S 

which means that 


-an-^Q 


nOO pOO 


n \ n 


/O JO 


exp -E^i ^ n (3.40) 

\ 7=1 / lc=l l^i<j^n 7=1 


Therefore 


l(s) :==5f(l) •X(”)(a,^) 

f-ocn+^Q-l 

Letting f = 1 gives the conclusion: 


(3.41) 


J(f) = ^-'(^(1)) = ^-'(1) = ■XW(«,/5). (3.42) 

T (an+ ^( 2 )) 




(3.43) 


In order to calculate the integral (a, j6), it suffices to calculate the integral (a, jS), which is 
derived in Corollary 16.41 via Selberg's integral (See Appendix): 


n r(l+;f)T(a+(;-l)f) 


/=i 


r(i + f 


(3.44) 


Matrix integrals, especially over unitary groups, are very important. We recently present 
some results of this aspect IIT^ . Generally speaking, the computation of matrix integrals is very 
difficult from the first principle. Thus frequently we need to perform variable substitution in 
computing integrals. The first step in substitution is to compute Jacobians of this transformation, 
this is what we present. We also apply the matrix integrals over unitary groups to a problem 
1201 l2n in quantum information theory. 


4 Applications 

4.1 Hilbert-Schmidt volume of the set of mixed quantum states 

The present subsection is directly written based on l24l l25l . The results were already obtained. 
We just here add some interpretation, from my angle, about them since the details concerning 
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computation therein are almost ignored. Any unitary matrix may be considered as an element 
of the Hilbert-Schmidt space of operators with the scalar product (U, V')hs = Tr(u*Vy This 
suggests the following definition of an invariant metric of the unitary group U{n)\ denote dG := 
U*dU, then 


ds^ := (dG, dG = - Tr ( dG^ ), 


implying 


ds2= E 

hM 


n ^2 ^ 

dG 


= E 

7=1 




2 ” _ 2 

+ 2E 




dGij 


Since dG* = —dG, it follows thaf 


ds^ = E 

7=1 


dG,, 


' + (d(Re(Gy)))' + 2^ (d(Im(Gy)))' 

i<j i<j 


(4.1) 

(4.2) 


(4.3) 


This indicates that the Hilbert-Schmidt volume element is given by 


that is. 


n(n-l) 


dv = 2^ "'nd (^(Gj,)) X nd{Re(G„))d(Im(G„)) = 2=^[dG], 

7=1 i<i 


voIhs (Z^(ii)) := [ dv = 2"'”2'’/' [dG] 

Ju{n) Ju{n) 


»(»-!) 

= 2 2 X 


U[n) 


n(n+l) 

Tn—r- 


1!2!--- (n-1)! 


n(n+l) 

_ (271) 2 

1!2!--- (m-1)!' 

Finally we have obtained the Hilbert-Schmidt volume of unitary group: 


(4.4) 


(4.5) 

(4.6) 

(4.7) 


voIhs = 2^vol {U{n)) = 


(4.8) 


We compute the volume of the convex (n^ — 1)-dimensional set D (C”) of density matrices of 
size n with respect to the Hilbert-Schmidt measure. 

The set of mixed quantum states D (C”) consists of Hermitian, positive matrices of size n, 
normalized by the trace condition 

D (C") = {p : C" ^ C'^l^ = p,p ^ 0,Tr (p) = 1}. (4.9) 
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It is a compact convex set of dimensionality (n^ ~ !)• Any density matrix may be diagonalized 
by a unitary rotation, p = UAU*, where U E U{n) and A = diag(Ai,..., A„) for Ay G R+. 
Since Tr {p) = 1, it follows that Ay = 1, so the spectra space is isomorphic with a (n — 1)- 
dimensional probability simplex A„_i {pG]R+:L;LiP; = i}. 

Let B be a diagonal unitary matrix. Since p = UBAB*U*, in the generic case of a non¬ 
degenerate spectrum (i.e. with distinct non-negative eigenvalues), the unitary matrix ti is deter¬ 
mined up to n arbitrary phases entering B. On the other hand, the matrix A is defined up to a 
permutation of its entries. The form of the set of all such permutations depends on the character 
of the degeneracy of the spectrum of p. 

Representation p = UBAB*U* makes the description of some topological properties of the 
(n^ — 1)-dimensional space D (C") easier. Identifying points in A„_i which have the same com¬ 
ponents (but ordered in a different way), we obtain an asymmetric simplex A„_i. Equivalently, 
one can divide A,j_i into n\ identical simplexes and take any one of them. The asymmetric 
simplex A„_i can be decomposed in the following natural way: 

A„-i = U .4, (4.10) 

ci\ -\ - \-dj^=n 

where k = 1 ,...,n denotes the number of different coordinates of a given point of A„_i, di 
the number of occurrences of the largest coordinate, ^2 the number of occurrences of the sec¬ 
ond largest etc. Observe that is homeomorphic with the set Gjt, where Gi is a single 

point, G 2 is a half-closed interval, G 3 an open triangle with one edge but without corners and, 
generally, Gk is an {k — 1 )-dimensional simplex with one {k — 2 )-dimensional hyperface without 
boundary (the latter is homeomorphic with an {k — 2)-dimensional open simplex). There are n 
ordered eigenvalues: Ai ^ A 2 ^ ^ A„, and n — 1 independent relation operators "larger(>) or 

equal(=)", which makes altogether 2”^^ different possibilities. Thus, A„_i consists of 2”^^ parts, 
out of which (”Ti) parts are homeomorphic with Gm, when m ranges from 1 to n. 

Let us denote the part of the space D(C”) related to the spectrum in ^ 4,,..,4 (k different 
eigenvalues; the largest eigenvalue has di multiplicity, the second largest ^2 etc) by A 

mixed state p with this kind of the spectrum remains invariant under arbitrary unitary rotations 
performed in each of the dy-dimensional subspaces of degeneracy. Therefore the unitary matrix 
B has a block diagonal structure with k blocks of size equal to di,... ,d]^ and 

D 4 . 4 ~ [ZT(n)/(W(di) X • • • X W(4))] X G,, (4.11) 

where d^ + ■ ■ ■ + dj^ = n and dj > 0 for j = 1,... ,k. Thus D (C”) has the structure 

D(C'^)~ U D 4 . 4 ~ U [W(n)/(W(di)x...xZT(4))]xG,, (4.12) 

di-\ - \-dj^~n rfpH- \-d)(=n 
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where the sum ranges over all partitions (di, • ■ • ,d^) h n of n. The group of rotation matrices B 
equivalent toU{di) x ■ • ■ xU (djt) is called the stability group of U (n). 

Note also that the part of represents a generic, non-degenerate spectrum. In this case 

all elements of the spectrum of p are different and the stability group is equivalent to an n-torus 

Di.1 ~ [U{n)/{V({iy^] X Gn. (4.13) 

The above representation of generic states enables us to define a product measure in the space 
D (C”) of mixed quantum states. To this end, one can take the uniform (Haar) measure on ZT(n) 
and a certain measure on the simplex A„_i. 

The other 2"^^ — 1 parts of D (C") represent various kinds of degeneracy and have measure 
zero. The number of non-homeomorphic parts is equal to the number P{n) of different represen¬ 
tations of the number n as the sum of positive natural numbers. Thus P{n) gives the number of 
different topological structures present in the space D (C”). 

To specify uniquely the unitary matrix of eigenvectors U, it is thus sufficient to select a point 
on the coset space 

■.= U{n)/U{l)^'^, 

called the complex flag manifold. The volume of this complex flag manifold is: 


voIhs (fIc 


voIhs {P({n)) ^ (2/1) 

vo1hs(Z^(1 ))” l!2 !..-(n-l)!- 


(4.14) 


The generic density matrix is thus determined by (n — 1) parameters determining eigenvalues 
and (n^ — n) parameters related to eigenvectors, which sum up to the dimensionality {n^ — 1) 
of D (C”). Although for degenerate spectra the dimension of the flag manifold decreases, these 
cases of measure zero do not influence the estimation of the volume of the entire set of density 
matrices. In this subsection, we shall use the Hilbert-Schmidt metric. The infinitesimal distance 
takes a particularly simple form 

dsRS = II dp IIhs = (dp. d^HS (4.15) 

valid for any dimension n. Making use of the diagonal form p = UAW, we may write 

dp = U (dA + [LfdLf, A]) U* (4.16) 
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Thus the infinitesimal distance can be rewritten as 


dsRS = 

i<i 

= EdA2 + 2f(A,--Ay) 






U*dU 


dG 


(4.17) 

(4.18) 


where dG = U*dU. Apparently, dAy = 0 since X!i/Li Ay = 1. Thus 




n—1 


dsRS = E dA/mydAy + 2E(A/ - Ay) 
'4=1 '</ 


dG, 


(4.19) 


The corresponding volume element gains a factor det(M), where M = [my] is the metric in the 
(n^ — n)-dimensional simplex A„_i of eigenvalues. Note that 


M — 11« + 


1 1 
1 1 


1 

1 


1 1 ■■■ 1 


Therefore the Hilbert-Schmidt volume element is given by 


dVHS = V^U dAy Yl{\i - Ay)2 


/=! i<j 


Ylld (Re(Gy)) d (im(Gy)) 


Then 



= ^/^2^vol(D(C'’)). 


(4.20) 


That is, respect to Hilbert-Schmidt measure, the volume of fhe set of mixed quanfum sfates is 

voIhs (D (C'^)) = (D (C’^)), 


i.e. 


voIhs (D (C")) = y/n{2Tz) 


^ r(i)r(2)---r(n) 

T(n 2 ) 


(4.21) 
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We see from the above discussion that the obtained formula of volume depends the used 
measure. If we used the Hilbert-Schmidt measure, then we get the Hilbert-Schmidt volume of 
the set of quantum states ll^ : if we used the Bures measure, then we get the Bures volume of 
the set of quantum states MT^ . 

A special important problem is to compute the volume of the set of all separable quantum 
states, along this line, some investigation on this topic had already been made Il2^l^ . There are 
some interesting topics for computing volumes of the set of some kinds of states, for instance, 
Milz also considered the volumes of conditioned bipartite state spaces [131, Link gave the ge¬ 
ometry of Gaussian quantum states IZl as well. We can also propose some problems like this. 
Consider the following set of all states being of the form: 

C := ||0i2 e D (C^ 0C2) : Tri(pi2) = ^l2 = Tr2(pi2)| . 

Parafharathy characterized the extremal points of this set 1141 . He obtained that all the extremal 
points of this convex set is maximal entangled states. That is, 

|0)|i/^o) + |l)|i/;i) 

It remains open to compute the volume of this convex set C 2 I 2 ) • 

4.2 Area of the boundary of the set of mixed states 

The boundary of fhe set of mixed states is far from being trivial. Formally it may be written as a 
solution of the equation 

det(| 0 ) = 0 

which contains all matrices of a lower rank. The boundary 3D (C'') contains orbits of different 
dimensionality generated by spectra of different rank and degeneracy. Fortunately all of them 
are of measure zero besides the generic orbits created by unitary rotations of diagonal matrices 
with all eigenvalues different and one of them equal to zero; 


A = {0,A2 < ■ ■ ■ < A„} . 


Such spectra form the (n — 2)-dimensional simplex ts.n- 2 , which contains (n — 1)! the Weyl 
chambers—this is the number of possible permutations of elements of A which all belong to 
the same unitary orbits. 
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Hence the hyper-area of the boundary may be computed in a way analogous to (|4.21ll : 


/ rank(X)=n—1 


[dX] = 


i.e. 


/ 0 <A 2 <---<A„ 


'EA-i) n |A/-A,fn(ydA,)x/' |dG,| 

/ 2=;<;<n ;=2 4Wi(«) 


n—1 

poo poo / 


(n - 1)! Jo 


MEA,-1 n |A,-A;|"n(AfdA,)x/ IdG 

\j=2 J 2=i<i^n j=2 JUi(n] 


n{n—1) 


1 r(i)---r(n)r(i)---r(n-M) 


Tin) 


r(n2-l) 


r(i)...r(n) 


n(n-i)rfH---rfn + H 


TT^ E/ (d fc^ 


r(n)r(n2-l) 


vol[^g = y/n — 12 2 vol*^” (t) = \/n — 1(2 


In)^ 


”("-!) r(i) • ■ ■ r(n -i- 1 ) 
^ r(n)r(n2-l) ■ 


In an analogous way, we may find the volume of edges, formed by the unitary orbits of the 
vector of eigenvalues with two zeros. More generally, states of rank N — n are unitarily similar 
to diagonal matrices with n eigenvalues vanishing. 


A - {A^ - ■ ■ • - A;JJ - 0,Affi^l < • • • < A,;} . 


These edges of order m are — m^ — 1 dimensional, since the dimension of the set of such spectra 
is n — m — 1, while the orbits have the structure of U(n) / {U(m) x U(1)"^’”) and dimensionality 
— {n — m). We obtain the volume of the hyperedges 


vol 


(n—m) 

HS 


1 ™Ihs (Hc’) 
(n - my. c(«".A) 


4.3 Volume of a metric ball in unitary group 


Consider a metric ball around the identity in the n-dimensional unitary group ZT(n) with 
Euclidean distance e, 


Be := {li G W(n) 


U-l„ 



(4.22) 
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where ||*||p is the p-norm for p = 2. We consider the invariant Haar-measure p, a unform 
distribution defined over lA{n). Denote the eigenvalues of 11 by The joint density of the 

angles 6j is given by 


p(0i,.. .,0„) = 


_ - _ n 

{lTzYn\ j 




(4.23) 


where 9j G [—tt, tt], / = 1,..., n. In what follows, we check the correctness of the integral formula: 

f p{9)d9 = 1. 


Indeed, set J{9) = Yli<ii<j^n 




and so 


i<j i<j 

= (signT)(^i • • ■ 

i<i 

where t = (n • ■ ■ 21), i.e. T(y) = n + 1 — / or t is written as 


(4.24) 

(4.25) 


T := 


1 2 
n n — 1 


n 

1 


n{n—l) 


Note that signr = (—1) 2 . We see that the integral is the constant term in 


C„(signT)(^i ■ • • (4-26) 

i<i 

Thus our task is to identify the constant term in this Laurent polynomial. To work on the last 
factor, we recognize 


V{0=Yli^i-Q) 


(4.27) 


as a Vandermonde determinant; hence 

ViO = x: (signu)C^(')-' • • • (4.28) 

c^Sti 


Hence 


n(?, - = V(Cf = E (sign^)(sign 

i<i a.TceSn 


Let us first identify the constant term in 


/(0) = (signT)(^i---^„)“Mv(O'. 


(4.29) 


(4.30) 
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We see this constant term is equal to 


{■In pin 



= (signT 


pin pin / 


■cr(l) + 7r(l) —n—1 ^cr{n)+n{n)~n—l 


de 




= (signT) Y, (sign (sign tt) — / Ci 

<T,neSn 

= (signT) ^ (sign tr) (sign tt) = (sign t) Y (sign cr) (sign tt) 

{a,n)e.Sny-S„->^j,a(j)+n(j)=n+l (o-,n)eS„xS„:n=ro- 

Note that the sum is over all (cr, n) E S„ x S„ such that cr{j) + 7i{j) = n + 1 for each j E {1,... ,n}. 
In other words, we get Ti{j) = n + 1 — cr{j) = T{a{j)) for all j E {1,. , n}, i.e. n = ra. Thus the 
sum is equal to 

(signr) E (signer) (sign Tu) = n\, 


aeSn 


which gives rise to 


1 /'27j: pin I 

Jl^Io '"lo = 


Now the condition on the distance measure 


U-ln 


^ e is equivalent to 


E 


_ J 


^ e^. 


Using Euler's formula = cos 9 + ^/—l sin 6 and the fact that (cos 0 — 1 )^ + sin^ 0 = 4 sin^ (^), 
we get 


U-% 


^ e 


M 


Thus the (normalized) volume of the metric ball equals the following: 

vol(Be) := ^{Be) = [ djiiU). 

JBe 

By specfral decomposition of unitary matrix, we have 


(4.31) 


(4.32) 


U = VDV*, D = e^®, 0 = 


01 0 
0 02 

0 0 


0 

0 

0 

9„ 


, 0y G [-7r,7r](; = 1,.. .,n). 


d0„ 
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Hence 


V* ■dU-V = dD 


V*dV,D 


implying that 

V* ■U*dU-V = D* (^dD+ V*dV,D 
Let dG = U*dU,dGi = V*dV and dX = V* ■ dG ■ V. Thus 

[dX] = [dG] 

because of V G U(n)/U(l)^". We also have 

dX = D* • (^dD + dGi, D ) =D* -V* -dU-V, 

it follows that 


(4.33) 


[dX] = [dU]. 


Apparently, 


Be = {yOW G U{n) : D G U{n), V G U{n), 


D-l„ 


^ e\ 


But 


[du]= n 

10<j^n 




[dD][dGi], 


therefore 


(4.34) 


(4.35) 


(4.36) 


[dG]= n 




[dD][dGi], 


(4.37) 


together with the facts that the region in which (0i, ...,6n) lies is symmetric and V E U(n) /lA 
implying 

2 


L n! / ■ ■ ■ / n 




,G^6i _ J^Bj 


[dD] X / [dGij. 
Ju{n) 


That is. 


d^{U) = 


(iTiYnl J "J 

7 . 9 , . ,_s ^ c 


n 

E"=isW(e^/2)s:f 


^G^Bi _ ^G^Bj 


[dD], 


where 


iu{n)\dG\ 7W(n) JU^{n) 
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From this, we get 


vol 


(Fe) f j P(^l/ ■ ■ ■ / ^n) 


d0,-. 


(4.38) 


where e G [0,2-^?!]. For the maximal distance e = 2^n, the restriction sin^(0;/2) ^ ^ 


becomes irrelevant and vol ( ^ 2 ^] = 1- 


We start by rewriting the n-dimensional integral (I4.38II . with the help of a Dirac delta function. 


as 


vol 

We know that 


/ oo 

5{t — a)f{t)dt = f{a). 

-OO 


p( 01 ,. . ., 0 „)]^d 0 ydf. 

7=1 


If /(f) = then 

/•P /-^ 

/ <5(f — fl)df = / <5(1 — fl)l{t:fl:^t<^}df = 

J DL J — OO 

By using the Fourier representation of Dirac Delta function 

1 /»00 _ 

we get 


l‘Ht-a)d, = i 


e^d-“)Mf ds 


1 

271 


-00 \ jO 

00 


-00 se 


V —IflS 


ds. 


Let a = W/i sin^ — 'i ~ I^f=i 2 Indeed, 


-£sin2(0y/2) = -- + ^(l-sin2(0y/2)) =-- + ^cos2(0y/2) 

7=1 7=1 7=1 

7=1 7=1 


Therefore 


= J?/ 




(4.39) 


7=1 


27r J —OQ ^ 2 ^ 


,_ ^ COS 0 J 


1 

In 




,_ COS 0 ,- 


-00 ^ 2 ^ 




ds 


V7=l 
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Inserting this formula into (I4.39I) and performing the integration over t first, we have 


1 .00 

vol{Be) = 7^ - 

2.7T J— 00 S€y 


1 

271 , 


-00 5 ^v ^2 


— 


n _ cos 6; 

/ , P{Ql .0n)n^'^'^d0y 

'[-7r,7r]» 


D„(s)ds, 


where 


Since 


Dn{s) = [ p{9l,- ■ ■,0n)Y\ 

4[-7T,7t]« 


,_ COS 6. 


\/^S- 


-dOi. 


gvA^e,' _ g\/~d-Bj _ ^g|. ^gvrT(! 1)6;;^ 

is a product of two Vandermonde determinants where i,k = 1,... ,n. The following fact will be 
used. 

Proposition 4.1 (Andreief's identity). For two nxn matrices M(x) and N(x), defined by the following: 



Mi(xi) 

Mi(x2 ) • 

A4i(x,j) 


Ni(xi) 

Ni(X2 ) • 

• Ni(x„) 

M(x) = 

M2 (Xi) 

M2 (x2 ) • 

■ MliXn) 

,N{x) = 

N2 (Xi) 

N2 (x2 ) • 

■ N2 (x,0 


M„(xi) 

M„(X2) • 

hAni^Xn) 


N„(xi) 

N„(x2 ) • 

Ain ( Xji ) 


and a function iv(-) such that the integral 

nb 

/ Mi{x)Nj{x)w{x)dx 
J a 

exists, then the following multiple integral can be evaluated as 

[■ ■ ■ I det(M(x)) det(N(x)) w(xy)dxj = det r/" Mi{t)Nj{t)w{t)dt\ , (4.40) 

^a,b j—\ \J a J 

where ts.a,h '■= {x = (xi,..., x„) G R" : fc ^ xi ^ X 2 ^ ^ x„ ^ a}. 

By invoking this identity, we know that 

D„(s) = ^-^[...[ _^det(c^('-i)''^)det(e-'^('-i)''^^^ 


.J [—n,n]’' 


7=1 


{2n)"nl J 

= 1 del ( /" . 

(27r)” V 7 -^ J 

In what follows, we need the Bessel function of fhe first kind. The definition is 
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g/^(n0-xsm0)^g _ 

M 


(-ly 

r(; + n + 1);! 


(.)2,« 


We can list some properties of this Bessel function: 

(i) g\/^ine+xsine)^g _ ^y/^{-n0-xsme)^0 

(ii) J„(x) = ^ ^ 

(iii) /_„(x) = ^ ^ ^ ^^(ne+xsme)^^ 

We claim that 


r ^V^{n0+xcos9)^0 ^ 2ne^'-^Ux). 

(4.41) 

J —71 



Indeed, since 

r ^^{ne+xcose)^Q ^ r e^{ne+xsme)^Q ^ g-v^'r . 27 zJ^„{x) 

J— 71 J— 71 

and ]~n{x) = (—l)”/„(x), the desired claim is obtained. 

Thus 

J — 7T J— 71 


= 2ne 




"Jh ( 2 ) • 


We now have 


D„(s) - det ( 2ne^'^^Ji 

>n > (2Ti)n ' 


- 1 C 2 


{2tz)^ 


de.(,,.,©). 
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By the definition of Jn{x), we have { — and furthermore. 


Dn{-s) = D„{s). 


Therefore 


Vol(Be) = 




-i^s 


'2.7Z J—oo S€^ 


D„{s)ds 


TC Jo s 

where we used Euler's formula. In fact, we have 


1 .a.sin(f) +sin((^-fjs, 

-I - 




se^ 


cos 


(fs) -cos(( 


T~i)y ^-sin(fs)-sin 
s s 

Finally, we get the volume of a metric ball in unitary group U{n): 


ds. 


T-Il« 



5 Appendix I 

5.1 Matrices with simple eigenvalues form open dense sets of full measure 

The present subsection is written based on the Book by Deift and Gioev [H. 

Let Jfn{C) be the set of all n x n Hermitian matrices with simple spectrum, i.e. the multiplic¬ 
ity of each eigenvalue are just one. Next we show that is an open and dense set of full 

measure (i.e. the Lebesgue measure of the complement is vanished) in H (C”), the set of all n x n 
Hermitian matrices. 

Assume that M is an arbitrary Hermitian matrix in with simple spectrum < 

■ ■ ■ < jin, then by standard perturbation theory, all matrices in a neighborhood of M have simple 
spectrum. Moreover if H G H (C”) with eigenvalues {hj n}, then by spectral theorem. 
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H = UAU* for some unitary U and A = diag(/zi ,... ,h„). Now we can always find Cj arbitrarily 
small for all j so that hj + ej are distinct for all j. Thus 

He := lJdiag(/ri + €i,... ,hn + €n)U* 

is a Hermitian matrix with simple spectrum, arbitrarily close to H. The above two facts show 
that is open and dense. In order to show that is of full measure, consider the 

discriminant: 

A(A):= n 

By the fundamental theorem of symmetric functions, A is a polynomial function of the ele¬ 
mentary symmetric functions of the Ay's, and hence a polynomial function of the entries 

Hu,..H„„, Re(H,y), Im(H;y) {i < j) of H. Now if (C) were not of full measure, then A would 

2 2 

vanish on a set of positive measure in K” . It follows that A = 0 because A is polynomial in K” . 

Let H = diag(l,2 ,... ,n) is a Hermitian matrix with distinct spectrum and A(H) 7 ^ 0. This 
gives a contradiction and so is of full measure in H (C”). 

Similarly, the set A^, 7 (R) of all n x n real symmetric matrices with simple spectrum is an open 
and dense set of full measure in the set S (R”) of all real symmetric matrices. 

From the above discussion, we conclude that the set of all density matrices with distinct 
positive eigenvalues is an open and dense set of full measure in the set of all density matrices. 

We begin by considering real symmetric matrices of size n, which is the simplest case. Because 
A^„(R) is of full measure, it is sufficient for the purpose of integration to restrict our attention 
to matrices M G A^„(R). Let }ii < ■ ■ ■ < }in denote the eigenvalues of M. By spectral theorem, 
M = UAir, where A = diag(/^i,.. Clearly the columns of U are defined only up to 

multiplication by ±1, and so the map M i-G (A, LI): 

^„(R) 9 M ^ (A, U) G R” X 0{n) 

is not well-defined, where 

• • ■ 'V-n) G R” : < • • • < , 0{n) = n x n orthogonal group. 

We consider instead the map 

(pi : =r„(R) 9 M ^ (A, Q) G R” X {0{n)/Ki) (5.1) 

where Ki is the closed subgroup of 0{n) containing 2" elements of the form diag(±l,..., ± 1 ), 
0(n)/K\ is the homogeneous manifold obtained by factoring 0{n) by Ki, and U = UKi is the 
coset containing U. The map (pi is now clearly well-defined. 
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The differentiable structure on 0{n')/K\ is described in the following general result about 
homogeneous manifolds: 

Proposition 5.1. Let K be a closed subgroup of a Lie group G and let G/K be the set {gK : ^ G G} 0 / 
left cosets module K. Let tz : G ^ G/K denote the natural projection n{g) = gK. Then G/K has a 
unique manifold structure such that n is G°° and there exist local smooth sections of G/K in G, i.e., if 
gK G G/K, there is a neighborhood W of gK and a C“ map t : W —)■ G such that tt o t = id. 


In other words, for each gK, it is possible to choose a g' G gK C G such that the map 
gK g' = rjgK)) is locally defined and smooth. 


For example, if G = IR^, the multiplication group of nonzero real numbers, and if K is the 
subgroup {±1}, then G/K = {x > 0} and n{a) = |fl| = {±fl}. If gK = {±fl} for some fl > 0, 
then T({±fl}) = a. Also, if G = 0(2) and Ki consists of four elements of the form 


±1 0 
0 ±1 




Mil Ui2 

U21 U22 


G 0(2), 


then 


gKi = 



Mil U12 


1 

CM 

S 

1 

1_ 


1 

CM 

1 

s 

1_ 



U21 U22 

r 

— M21 U22 

t 

M21 —U22 

r 


— Mil —Ul2 

— U2I —U22 


Since m|^ 


^21 = 1 = ^12 + 1 ^ 22 , each column of g contains at least one nonzero number: For 
example, suppose M 21 and U 22 are nonzero; then there exists a unique g' G gK such that the 
elements in the second row of g' are positive. The same is true for all g"K close to gK. Then 
g"K I— )• g” is the desired (local map) t. We can generalize the above construction on 0{nj/K\. 
We will prove the following result: 


Proposition 5.2. (pi is a diffeomorphism from onto x {0{n)/Ki). 


Proof. Note that (here and below we always speak of the real dimensions) 

dim(R!( X (G>(n)/Ki)) (5.2) 

as it should. 

Define the map ^1 : R" x {0{n)/Ki) — t ynjlR) as follows: If (A, U) lies in R" x {0{n)/Ki) — > 
,i^„(R), then 

4>i{A,U) = UAW (5.3) 

where U is any matrix in the coset U. If U' is another such matrix, then U' = Uh for some h ^ Ki 
and so 

U'AjU'y = UhAh'^W = UAW. 
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Hence is well-defined. We will show that 

(/)1 o = id] 


(5.4) 


R1x{0{n)/Ki) 


and 


0^1 = id^^(K). (5.5) 

Indeed, 

(pi (^Am)) = (pi{UAW), ueu (5.6) 

= (A,LfXi) = (A,a) . (5.7) 

Conversely, if M = UAW G then 

^i{(pi{M)) = (^A,U= UKi'j = UAW = M. (5.8) 

This proves Eq. (15.41) and Eq. (I5.5E In order to prove that (pi is a diffeomorphism, it suffices to 
show that (pi and ^i are smooth. 

The smoothness of (pi follows from perturbation theory: Eix Mq = UoAqUq G =5^„(R), where 
Ao = diag(/ioi/• • • /l^On) for }ioi < < jion, and Uq G 0(n). Then for M near Mq, M G 

^„(R), the eigenvalues of M, f(i(M) < ■ ■ ■ < are smooth functions of M. Moreover, the 

associated eigenvectors Uj(M), 

Muj{M) = \j{M)uj{M), l^i ^n, 


can be chosen orthogonal 


{ui{M),Uj{M)) = 5ij, (1 ^ i,i ^ n). 


and smooth in M. Indeed, for any j with 1 ^ ^ n, let Pj be the orthogonal projection 

Iff 




s-M 


-ds 


(5.9) 


where Pj is a small circle of radius e around }iQj, | }ioi ~ Eo/1 > £ for i 7 ^ j. Then, where wy(Mo) is 
the ;-th column of Mq, 


Uj{M) = 


Pj{M)uj{Mo) 


{ujiMo),Pj{M)ujiMo)) 


1 ^ 7 ^ n. 


is the desired eigenvector of M. 
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Set U(M) = [ui(M),.. .,u„(M)] E 0{n). Then M i—)■ (A(M),U(M)) is smooth and hence 

(pi(M) = (A(M),fJ(M)) = ti{U{M)) 

is smooth, as claimed. 

Finally, we show that is smooth. Fix ^Aq, G R” x {0{n)/Ki) and let r be the lifting 
map from some neighborhood W of Uq to 0{n). Now for all U E W, t{U) E U by Proposition l5.ll 
Hence 

(^A,u) - T{U)AT{Uy (5.10) 

from which it is clear that is smooth near ^Aq, Lfo^, and hence everywhere on R” x (0(n)/Xi). 
This completes the proof. □ 

Now consider J^(C). The calculations are similar to =5^„(R). Define the map 

(p2 : 3 (A,Lf) e R^ X {U{n)/K2). (5.11) 

Here K 2 is the closed subgroup oiU{n) given by T x • • ■ x T = {diag ..., G R}, 

A = diag(Ai,..., A„), Ai < • • • < A„ as before, and 

M = UAU*, U = 7z(U), 

where n : lAin) -3 U{n)/K 2 is the natural projection U i-G 7z(U) = LIK 2 of the unitary 
group U{n) onto the homogeneous manifold U{n)/K 2 . Let t : U{n)IK 2 —t U{n) denote the 
(locally defined) lifting map as above. As before, (p 2 is a well-defined map that is one-to- 
one from J^n{C) to U{n) -3 U{n)/K 2 . Note that because dim(ZL(n)) = ,diTa{U{n)/K 2 ) = 

dim.{lA{n)) — dim(X 2 ) = — n, as it should. 

As before, the proof of the following result is similar to that of Proposition 15.21 

Proposition 5.3. (p 2 is a diffeomorphism from Jf^(C) onto R" x {U{n)/K 2 ). 

5.2 Results related to orthogonal groups 

Proposition 5.4. Let Y,X,T E R”^", ivhere Y and T are nonsingular, X is skew symmetric and T is 
lower triangular of independent real entries yifs, Xjj's and tifs respectively. Let tjj > 0,j = l,...,n — 
1,-00 < t„„ < 00 ; —00 < tjk < co,j > y —00 < Xjk < CO,j < k or —00 < tjk < co,j ^ k and 
the first row entries, except the first one, of {1„ -|- X)^^ are negative. Then the unique representation 
Y = T [2{ln + X)“1 - In] = T(ln - X) (1„ + X)-i implies that 

[dY] = (det(ll„ + X))“(”“^) [dr][dX]. (5.12) 
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Proof. Take the differentials to get 

dY = dT-+ +T- (-2(l„ + X)-i-dX-(l„ + X)-i) 

= dT ■ (2(1„ + X)-i - In) + T ■ + Z) ■ dX • (1„ + Z)^ 

= dT-Z-iT(l„ + Z)-dX-(1„ + Z). 

where Z = 2(1,, + X)“^ — 1„. Then we have 

P-^ • dY • Z^ = P-^ ■ dP - i(l„ + Z) ■ dX • {In + Z^). 

The Jacobian of the transformation of T and X going to Y is equal to the Jacobian of dT, dX going 
to dY. Now treat dT, dX, dY as variables and everything else as constants. Let 

dli = T-i • dY ■ Z^ dV = r-i • dT, dW = (1„ + Z) • dX • (1„ + Z^). 

Thus 


[dU] = det(T)-Met(Z^)'’[dY] = det(T)-'’[dY] 

^ [dY] =det(T)'’[dli] = (^n|f;7r) W- 

Note that det(Z^) = ±1 since Z is orthogonal. Since X is skew symmetric, one has 
[dW] = det(l„ + zy-^ldX] = det(l„ + X)-("-i)[dX]. 


One also has 


[dy] = 



[dT]. 


Now we see that 

dli = dY - -dW U = V --W. 

2 2 

Let U = [uij],V = [r;;y],W = {wif. Then since T is lower triangular tij — 0, z < i and thus V is 
lower triangular, and since X is skew symmetric Xjj = 0 for all j and x,y=_ 3 ;.. for i j and thus W 
is skew symmetric. Thus we have 

1 . . 

Uii = Vii,Uij = <] 

1 

Uij = Vij + -Wij,l> ]. 
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Take the w-variables in the order u,-,-, z = 1, ...,n;Wy,z < i;Uij,i > j and VuJ = 1,... ,n-,Vij,i < 
i',Vij,i > j. Then the Jacobian matrix is of the following form: 

1 „ 0 0 

0 0 

. 0 (2)1© 1© . 

and the determinant, in absolute value, is (I)”*"” That is, 

/I \ n(n-l)/2 

[du] = (^-j [dy][dw]. 

Now substitute for [dU], [dW] and [dV], respectively to obtain the result. □ 

Proposition 5.5. Let Y,X,D G he of independent real entries, where Y is symmetric with distinct 
and nonzero eigenvalues, X is skew symmetric with the entries of the first row of {1„ + X)~^, except 
the first entry, negative and D = diag(Ai,.. .,A„), with Ai > • • • > A„. Then, excluding the sign, 
Y = [2{1„ + X)-i - 1„] D [2{ln + X)-i - In] implies that 

[dY] = {det{ln + X))~‘'''~d ^n|A;-Ay|^ [dX][dD]. (5.13) 

Proof Let Z = 2(1„ + X)~^ — 11„. Take the differentials and reduce to get 

ZN dY ■ Z = -^(1„ + Z^) ■ dX ■ (1„ + Z) • D + dD + ■ (1„ + Z^) ■ dX • (1„ + Z). 

Put 

dU = Z^-dY-Z, dW = dD, d© = (1„ + Z^)dX(l„ + Z). 

Thus 

dU = -^dV • D + ■ dY + dW. 

But since dY is symmetric and Z is orthogonal, [dli] = det(Z)”+^[dY] = [dY], excluding the 
sign. Clearly [dW] = [dDj. Since X is skew symmetric we have 

dY = det(l„ + Z)”-i[dX] = 2 <’^-d det(l„ + X)-(”-i)[dX]. 


We see that 


duii = dwii, 

1 

duij = -{\i-\j)dvij,i < j. 
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Take the w-variables in the order ua, i = 1,... ,n; uij, i < j and the zv and z;-variables in the order 
Wii,i = 1,... ,n}Vjj,i < j. Then the matrix of partial derivatives is of the following form; 

’ 1 o' 

0 M J ' 

where M is a diagonal matrix with the diagonal elements — ^j),i < /• There are n{n — l)/2 
elements. Hence the determinant of the above matrix, in absolute value, is 2 n[n Ty 

That is. 


Hence 


[dii] - I a,- - Ay|^ [dV][dD]. 

[dY] = det{ln + [dX][dD]. 


□ 


Remark 5.6. When integrating over the skew symmetric matrix X using the transformation in 
Proposition 15.51 under the unique choice for Z = 2(1„ + X)“^ — 1„, observe that 

2 n{n-l)/2 f ^ X)-("-d[dX] = (5.14) 

Note that the Ay's are to be integrated out over oo > Ai > • ■ • > A„ and X over a unique choice 
of Z. 


5.3 Results related to unitary groups 

When X is skew hermitian, that is, X* = —X, it is not difficult to show that 1 ± X are both 
nonsingular and Z = 2(1 + X)^^ — 1 is unitary, that is, ZZ* = 1. This property will be made use 
of in the first result that will be discussed here. Also note that 

2(1 + X)“i - 1 = (1 + X)-i(l - X) = (1 - X)(l + X)-\ 

When X is skew hermitian and of functionally independent complex variables, then there 
are p + j-eal variables in X. Let T be a lower triangular matrix of functionally 

independent complex variables with the diagonal elements being real. Then there are real 
variables in T also. Thus combined, there are 2n^ real variables in T and X. It can be shown that 

Y = f (^2(1 + X)“i -l) 
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can produce a one-to-one transformation when the tjj's are real and positive. 

y = fz, 

where ZZ* = 1, T = tjj = tjj > 0,j = 1,... ,n. Then 

YY- = Tf =*(?, = f; lyitl". 

k=l 

Note that when tn is real and positive it is uniquely determined in terms of Y. Now consider the 
first row elements of TT* that is tntzi, ■ ■ ■, tnini- Hence f 2 i/ • • • / ini, that is, the first column 
of T is uniquely defermined in terms of Y. Now consider the second row of TT* and so on. 
Thus T is uniquely determined in terms of Y. But Z = T^^Y and hence Z, thereby X is uniquely 
determined in terms of Y with no additional restrictions imposed on the elements of Z. 

In this chapter the Jacobians will also be written ignoring the sign as in the previous chapters. 

Proposition 5.7. Let Y, X and T = [Ty^] ben xn matrices of functionally independent complex variables 
where Y is nonsingular, X is skew hermitian and T is loiver triangular with real and positive diagonal 
elements. Ignoring the sign, if 

Y = f (^2(X + l)“^-l) = f{l-X){l + X)-^, 

then 




det( (1 + X) (1 - X)) • [dX] [dT]. (5.15) 


Proof. Taking differentials in Y = T ^2(1 -|- X) ^ — 1^, one has 

dY = f (^-2(1 + X)-i ■ dX ■ (1 + X)“i) + df ■ (^2(1 + X)-^ - 1) . 
Let 

Z = 2(1 + X)~^ - 1 ^ (1 + X)~^ = i(l + Z) 
and observe that ZZ* =1. Then 


dY • Z* = --(1 + Z) • dX ■ (1 + Z*) + T-^ ■ dT. 


( 5 . 16 ) 
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Let dW = (1 + Z) ■ dX • (1 + Z*). Then 

[div] = det((l + Z)(l + Z*)) [dX] =22”"- det((l + X)(l - X)) • [dX], 

Let dU = • df, then 


Let dV = T-i • dY • Z*, then 


Eq. H5.16II reduces to 


[diJ] - 1 1 ■ [dT]- 


d=i 


[dV] = det(rr) ■ [dY]. 


dV= --dW + dU. 


( 5 . 17 ) 


( 5 . 18 ) 


(5.19) 


( 5 . 20 ) 


Note that dW is skew hermitian. Denote 


y = [vjk] = 

^ = [^jk] = [^]k] + 

u = [Uik] = [u]l^] + V^[u]h 

where v^j^\ m = 1,2 are all real. Then from Eq. (15.201) , we see that V = — + U. 

Thus 


(m) 

^)k = 

+ u^^\j >k,m = 1,2, 

(m) 

V-, = 

!k 

~\^]k^'j < Km = 1,2, 

.d) = 

u^\ 

11 

11 ' 

(2) 

V-- = 

1 (2) 

11 

2 n 

The matrices of partial derivatives are 

the following: 


9 

(vf,vf-,v^^,vf{j < K)-,v^^,vf{j > k)') 



< k)-,uf^\u^l\i > k) 

In 0 0 0 0 0 



0 -ll„ 0 0 0 0 

0 0 -\\) 0 0 0 

0 0 0 -IUq 0 0 

0 0 2l(«) 0 1(«) 0 

0 0 0 ll(«) 0 1(«) 
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where 


( 1 ) 


All — 


A 22 = 


( 2 ) 


,J = l,...,n j = In, 




,J j = --1„ 


and 


^33 = 


A 44 — 


A55 = 


^66 = 


2^53 = 


dv) 


( 1 ) 


. V 


dv 


( 1 ) 

Jk 






dv 


( 2 ) 

jk 


dwf^ 


^ I = “2^© 


dv^l'^ 

dv^l^ 

-^yi >k \ = lyy. 


dll 


jk 


dv 




■■ j > k \ — -1(«), ^64 — 


( 2 ) 

. V 


a^f 


■i>k \ - -1(«). 


The determinants of An,..., Aee contribute towards the Jacobian and the product of the deter¬ 
minants, in absolute value, is = 2^"^. Without going through the above procedure 

one may note from (I5.16D that since dX has real variables, multiplication by ^ produces the 
factor 2^”^ in the Jacobian. From (I5.17I) , (I5.18I) , (I5.19D and (I5.20I) we have 


= n f ■ [d^] = n f ■ 2^”' • [dw] m 


v;=i 

2 n 

m 

v;=i 


v;=i 


2n \ . 22n^-n^ . det((l + X)(l - X)) 


v;=i 


T2;-l) 


2”' 


\/=l 

This completes the proof. 


det((l + X)(l-X)) -[dXJidT]. 


[dX][dr] 


□ 


Example 5.8. Let X G skew hermitian matrix of independent complex variables. Then 

show that 
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hdX] det((ll + X)(l-X)) 
Jx 


71 " 


2”(”-i)r„(n) 


Indeed, let Y = [yjk] he a n x n matrix of independent complex variables. Consider the 
integral 


/idf|.-'(»-) = n 


^ p-\-oo 


e dyjk = 7z" . 


(5.21) 


Now consider a transformation Y = T ^2(1 + X)“^ — 1^, where T = [Fyjt] is lower triangular with 
tjj's real and positive and no restrictions on X other than that it is skew hermitian. Then since 
Tr = Tr (^TT*^, from Proposition 15.71 and Eq. (15.2111 . we have 


2(n-y)+l 

a 


det((l + X)(l-X)) 


= ^[df]e'^''(^'^*) -2”^ X ^[dX^ det((l + X)(l - X)) 


Note that 


-Tr(TT*) _ 


W=1 


\ ;=i i>k 


But 


Hence 


/ + O 0 

' 

-00 


e dfyj; = n and 


r*+oo 


j2(n—j)+l —flj, l-r./' • I d\ 

t. e ndtjj = -T{n-] + l). 




Substituting this the result follows. 

Remark 5.9. When a skew hermitian matrix X is used to parameterize a unitary matrix such as 
Z in Proposition 15.7[ can we evaluate the Jacobian by direct integration? This will be examined 
here. Let 


.y„:= f[dX] det((l + X)(l-X)) 
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Partition 1 + X as follows: 


1 + X = 


1 + xn X 


12 


-X*2 1 + Xi 


where X 12 represents the first row of 1 + X excluding the first element 1 + x^, and I + Xi is 
obtained from 1 + X by deleting the first row and the first column. Note that 


det(l + X) = det(ll + Xi) (1 + xii + Xi 2 (l + Xi)- 1 X 1 * 2 ) ■ 


Similarly, 


1-X = 


1 — Xn —Xu 

XI 2 1 - Xi 


and 

det(l - X) = det(ll - Xi) (1 - In + Xi2(l + Xi)-1X1*2) ■ 

For fixed (1 + Xi) let Uu '■= Xi2(l + Xi)-i, then 

[dai2] = |det((l + Xi)(l - Xi)) 1 “' • [dXi2] 
and observing that X^ = — Xi we have 

Xi2(l + Xi)-lXi*2 = ai2(l - Xi)Lfi*2. 

Let Q be a unitary matrix such that 

Q*XiQ = diag(^Ai,..., v^A„_i) 

where Ai,..., A„_i are real. Let 

1^12 = UuQ = [v-l, . . . , Vn-l]. 

Then 

1/12(11 - Xi)Lfi *2 = (1 - y^Ai) |yi |2 + • • • + (1 - ^A„_i) | 5 -„_i |2 

and 

1 + Xii + Xi2(l + Xi)^^Xi*2 = a - 

where 

a = 1 + I 1^ + ■ ■ • + I u,7_i 1^ 

(' 2 .) I ~ 12 I ~ 12 

b = ~^ii +Ai|yi| + • ■ ■ + A„_i I u„_i I 
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observing that xn is purely imaginary, that is, Xu = where x'^ is real. Thus 

det((l + X)(l-X)) 


yields the factor 


[{a - y/^b){a +y/^b)] " = {a^ + b^y 


So 


= L L [ (det{l + Xy{a-V^b)det{l-Xy{a + ^b)) " [dXi] [dXi2]dTii 

JXi JXi2 Jxii V / 

= L L [ det((l + Xi)(ll-Xi)) + by-’ydXy[dXu]dxn. 

JXl JXl7 J-Cll V / 


Since 


[dyi 2 ] = [dUi 2 ] and [dUi 2 ] = det (1 + Xi) (1 - Xi) [dXi 2 ], 


it follows that 


-(n-l) 


+ '’[dXi][dyi2]dxii. 


[dXi 2 ]= det(^(l + Xi)(Il-Xi)) [dyi 2 ]. 

Based on this, we have 

= 4 4/ det((l + Xi)(l-Xi)) 

JXi JVi2 Jxn V / 

Then 

JVi2 JXy^ 

= (l + 5) dxP’ 


ni 
b^\~^ 


Consider the integral over , —oo < < +oo. Change to b and then to c = b/a. Then 


L 


1 + ^^ "dx^'^ - 
p) r ^ fl2) - 


':=fc 


= 2a /“(I + c2)-Mc = /( 2 )r^ 1) 

Jo r(n) 


by evaluating using a type-2 beta integral after transforming u = c^. Hence 


= k^n-1 L [dVi2\a 
J V]2 




= 


n—1 


-00 /* + oo 


-00 J —00 


I ~ |2 I ~ |2\ ~(2«-l) _ _ 

1+ ^1 + • ■ ■ + ) d^i • • ■ dVn-\. 
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For evaluating the integral use the polar coordinates. Let vj = where and 


( 2 ) 

Vj ^ are real. Let 


( 1 ) o 

Vj ' = Vj cos 9j, 

( 2 ) • a 
V- = Vj smdj, 


0 ^ r,- < oo, 0 ^ 6; ^ 2tt. 


Then denoting the multiple integral by we have 


-(2n-l) 


dz7i • ■ ■ dz; 


n-l = / ■■■/ + +''' + |^«-l| 

J—CO J—CO V 

= (2.r-r...r 

v/1-1=0 Jr„_i=0 


Evaluating this by a Dirichlet integral we have 

.n-l r(n) 


X „_1 = tt '” 


r(2n - 1) 


for n ^ 2. 


Hence for n ^ 2, 


iAj = -TyrTT- 

r(2n — l)r(n) 

By using the duplication formula for gamma functions 


r(2n - 1) = ^22"-2p ( M - - ) r(n). 


Hence 


^n — ^n—1 


n 


22n-2r(„) 


Repeating this process we have 


TT 


n‘ 


n—l 


71 


TT 


22"-2r(n) 22('^-l)-2r(p - 1) 22-2r(l) 2"("~l)f„(n)' 

This is what we obtained in Example 15.81 


dr„_i. 


Next we consider a representation of a hermitian positive definite matrix Y in terms of a skew 
hermitian mafrix X and a diagonal mafrix D such that 

Y = ( 2(1 + X)~i - 1 ) D (^2(1 + X)-^ - 1 )* 

where D = diag(Ai,..., Ap) with the Ay's real distinct and positive, and the first row elements of 
(1 + X)“^ real and of specified signs, which amounts to require 2(1 + X)^^ — 1 G W {n)/U ( 1 ) 

In this case it can be shown that the transformation is unique. Note that 

y = ZDZ* = AiZiZj* + ■ ■ • + ApZ„Z* 
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it indicates that 


{Y-\jl)Zj = 0,i = l,...,n 

where Zi ,... ,Zn are the columns of Z such that (Zy,Z^) = 3jk- Since M, ..., A„ are the eigen¬ 
values of Y, which are assumed to be real distinct and positive, D is uniquely determined in 
terms of Y. Note that Zy is an eigenvector corresponding to Ay such that (Zy, Zy) = 1,; = 1,..., n. 
Hence Zy is uniquely determined in terms of Y except for a multiple of ±1, If any par¬ 

ticular element of Zy is assumed to be real and positive, for example the first element, then Zy is 
uniquely determined. Thus if the first row elements of Z are real and of specified signs, which is 
equivalent to saying that the first row elements of (1 -|- are real and of specified signs, then 
the transformation is unique. 

Proposition 5.10. Let Y and Xben xn matrices of functionally independent complex variables such that 
Y is hermitian positive definite, X is skew hermitian and the first row elements of (1 -|- X)~^ are real and 
of specified signs. Let D = diag(Ai,.. .,A„), where the Ay's are real distinct and positive. Ignoring the 
sign, if 

Y = (2{X +1)“^ - l) D (2{X +1)“^ - l) *, 

then 


|dY] = 


n|A|.-A/| 

\j>k 


det((l + X)(l-X)) 


• [dX][dD]. 


Proof. Let Z = 2(l-|-X) ^ — 1, X* = —X. Taking the differentials in Y = ZD Z* we have 

dY = dZ • D • Z* + Z ■ dD ■ Z* + Z • D ■ dZ*. (5.22) 


But 


and 


dZ = -2(l + X)~^-dX-(l + X)“^ 
= -i(ll + Z)-dX-(1 + Z) 

dZ* = 2(l-X)-^-dX-(1-X)“i 
= ^(l + Z*)-dX-(1 + Z*). 
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From (15.2211 . one has 


Z*-dY-Z = -^(1 + Z*)-dX-(1 + Z)-D + dD + ^D-(1 + Z*)-dX-(1 + Z) 

observing that Z*Z = 1. Let 

dU = Z* • dY • Z ^ [dU] = [dY] since Z*Z = 1, 

dV = (1 + Z*)-dX-(1 + Z) = (1 + X*)-^-LdX-(l + X)-i ^ 

[dY] = 4"'■ det((l + X)(l-X)) • [dX] 


(5.23) 


(5.24) 


if there are free real variables in X. But in our case there are only — n real variables in X 
when X is uniquely chosen and hence 

[dY] = 4""-" det((ll + X)(ll-X)) "-[dX], (5.25) 

and 


dU = -idY-D + ^D-dY + dD. 


(5.26) 


From (15.261) and using the fact that dY is skew hermitian and dLf is hermitian we have 


dujj = dAj, dujf ^ = ±^(A^ - Aj)dv^^\i >k,m = 1,2. 


Thus the determinant of the Jacobian matrix, in absolute value, is 


Ki>k 


j>k 


That is. 


[dU] ■p]]|Afc-Ay|^- [dY][dD]. 

j>k 

Substituting for [dLf] and [dY] from (15.231) and (15.241) the result follows. 


□ 


Example 5.11. For X an n x n skew hermitian matrix with the first row elements of (1 + X) ^ 
real and of specified signs show that 


AdX] det((Il + X)(Il-X)) 
Jx 


T„{n) 
A ' 


where 




f [dD] 

n |A,-A,f 




-Tr(D) 
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with D = diag(Ai,..., A„), Ai > ■ ■ • > A„ > 0. Consider a nxn hermitian positive definite 
matrix Y of functionally independent complex variables. Let 


B = 


= [ [dY]e-^<^) = [ [dY] det(Y) " 

Jy=y*>o Jy>o 


= r„(n) = 7 r^r(n)r(n-l)---r(l) 

= 7r“^^^(n — l)!(n — 2)! ■ ■ ■ 1! 

evaluating the integral by using a complex matrix-variate gamma integral. Put 

Y = ZDZ*, Z = 2{1 + X)-^ -1 

as in Proposition ISdO] Then 

[dY] - ■ (^n|A,-A^f^ ■|det((l + X)(l-X))P'-[dX][dD] 


and 


B = 


= hdX] det((ll + X)(l-X)) 

Jx 

f [dD]2'’('’-i) ■ 

iAi>->A„>0 J 


-Tr(D) 


Hence the result. From (i) in Example 13.21 [ we see that 


A = ( - j (r.(n)) , 


which implies that, when X is taken over all skew hermitian under the restriction that the first 
row elements of (1 + X)~^ are real and of specified signs. 


lx 


[dX] det((l + X)(l-X)) 


-n /jixn(n-l) I 


Tn{n) 


Remark 5.12. In fact, we can derive the volume formula (13.361) from Proposition 15.101 The rea¬ 
soning is as follows: 


L [dr| = 

JY>0:Tr(Y)=l JAi>-->A„>0 j 

X f [dX] det((l + X)(ll-X)) 

Jx 




Ki>k 


T(l)---T(n)T(l)---T(n + l) /7r\«(«-i) 1 

n! "" f(W) "" V2 J f„(n) 

^M^r(l).-.T(n) 

= 7T 2 - 

r(n2) ■ 
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Example 5.13. Let X G be a hermitian matrix of independent complex variables. Show that 


lAdX]e 


Tr(XX*) ^ 


Indeed, X* = X implies that 


Tr(xx*j = x:4+2x:i%r- 

;=1 i<i 


Thus 


jAdXV 


-Tr(XX*) _ 



n n{n-l) n{n-l) n{n-l) n2 

= 712x2 2 XTIZ =2 2 77-2. 


Example 5.14. By using Example 15.131 or otherwise show that 


L 


oo>Ai>--->A„>—oo 


I- Ay M exp [ - Ay ] dAi • ■ ■ dA„ = 2 


n{n-l) n 
2 J12 


^j>k 


;=1 


n-1 

m 

/=1 


Indeed, in Example 15.131 letting X = UDU*, where D = diag(Ai,.. .,A„) and U G Ui{n), gives 
rise to 


[dX] = [dD][dGi], Ai > • • • > A„ 


which means that 


n(n-l) 

2 TT 2 = 


j^xy 


.Tr(XX*) 


- L 


oo>Ai>--->A„>—oo 


Yli^i - exp - ^ A^ [dD] X / [dGi]. 




7=1 


Ui(n) 


That is. 


L 


oo>Ai>--->A„>—oo 


n(A. - Ay)^) exp f- E A?1 fl^Ay = 2-^7r5 Hy!. 
\i<i ) \ M M 


Therefore 


n(A, - exp (- t A;^) fidAy = 2-^715 H;!. 

G<; / V 7=1 / 7=1 7=1 
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6 Appendix II: Selberg's integral 


This section is rewritten based on Mehta's book irTTH . The well-known Selberg's integral is calcu¬ 
lated, and some variants and consequences are obtained as well. 

Theorem 6.1 (Selberg's integral). For any positive integer N, let [dx] = dxi ■ • ■ dxjv. 


A(x) = A(xi,...,xn) ^ 


Xj), ifN>l, 
h ifN = l, 


( 6 . 1 ) 


and 


N 


0(x)^0(xi,...,Xn)= |A(. 


1 27 




Then 


N 1 

SNioc,^,7)= •■•/ 0 (x)[dx] = n 

Jo Jo T-n 


r(x + 7 ;)r(^ + 7 ;)t (7 +1 + 7 /) 
T(a + /3 + 7(N + y-l))T(l + 7)' 


and for 1 ^ K ^ N, 


( 6 . 2 ) 


(6.3) 


/,’■■■ I' (n^') =n /; ■ ■ ■ 

valid for integer N and complex a, f>, 7 with 

Re(a) > 0, Re(^) > 0, Re( 7 ) > - min 


(6.4) 


(6.5) 


Aomoto's proof. For brevity, let us write 


(/(xi,.. .,X]v)) := 


/[o,i]^/(^i' ■ ■ ■' ..., XN)dxi ■ ■ ■ dxN 

0(xi,..., XN)dxi • ■ ■ dxN 


Firstly, we note that 
d 


dxi 


(x^X 2 • • • Xj( 0 ) 


= {a + x- l)x" ^X 2 • • • X^O - (j 6 - 1 )^ 




^1 


“2 ^ 1 -^/ 


Indeed, 


^ {x{X2 ■ ■ ■ Xk<^) = ax\ ^X2 • • • XkO -6 X?X2 • • • Xk^O(xi, . . . , Xn), 


(6.6) 


(6.7) 

(6.8) 


(6.9) 
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where 


Since 



N 


lAwrn 


;=i 



1x^X2 ■■ ■ Xk^(xi, ..., = X2 • • • Xk^(1,X2,...,Xn) -0 = 0, 


it follows from integrating between 0 and 1 in (I6.7II . we get 


0 = (fl + a 


l)(xi ^X2 ---xk} - (^- 1 ) 


X^X2---Xk\ , ^ /x1X2---Xk\ 


Now for a = 1 and a = 2, we get 


0 

0 


Clearly 


= a{x2 • • • xk) - (j 6 - 1 ) 


XiX2---Xk\ , ,, ^ IxiX2---Xk 


= {tt + l){xiX 2 - ■■ Xk) - iji-l) 




/ 3:1X2 • • • \ _ / x^X2 ■ • • \ 

\ Xi-Xy / \ / 


/ Xi3:2 • • • _ ^1^2 ■■■Xk \ 

\ Xi- Xj Xi - Xj / 

{xiX2 ■ ■ ■ xk); 


interchanging Xi and Xj and observing the symmetry. 


and 


/ X 1 V 2 ■■■Xk \ ^ _ / XjX2 ■■■Xk \ 
\ Xi-Xj / \ Xi-Xj / 


0, if 2 ^ 7 ^ K, 

\{xi---Xy), iiK<i^N, 


/ xlx 2 ■■■xk \ 
\ Xi-Xj / 


xjx2 ■■■Xk 

Xl - Xj 


i(viX2 ■ • • Vx), if2^;^K, 
{xiX2 ■ ■ ■ Xk), if X < ; ^ N. 


Performing the difference: (I6.13D — (I6.14I) and using the above, we get: 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

(6.14) 

(6.15) 

(6.16) 

(6.17) 


(a + 1 + j{K - 1) + 27 (N - X) + j6 - l)(xiV 2 ■ ■ ■ Xk) = {x + j{N - K)){x 2 ■ ■ ■ Xk), (6.18) 


or repeating the process 


(xiV2 ■■■Xk) 


(a + 7(N-X)) 
a + ^ + 7(2N-X-l) 

A ft + 7(N-;) 
fJa + )6 + 7(2N-;-l)’ 


(xi• • • xk-i) = 


(6.19) 

( 6 . 20 ) 
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For K = N, (I6.4D can be written as 


Sn(^ + l/i6,7) 

SN(a, /3,7) 


K 

= {xiX2 ■ • • ^n)= n 
M 


a + 7 (N -;■) 
a + fi + j{2N - i - 1) 


or for a positive integer a, 

^ SN(ft/i6,7) S]v(2,/3,7) 

SNih^rl) ~ SN(a-l,j6,7) SN(l,i6,7) 

X tt-l + 7 (N-;) A 1 + 7(N-;) 

a - 1 + ^ + 7(2N - 1) ''' 1 + ^ + j{2N - j - 1) 

A_ (ft - 1 + 7(N -;•))■■■ (1 + 7(N - ;•)) _ 

(ft-l + /i + 7 ( 2 N- 7 - 1 )) ■■■(! +/i + 7 ( 2 N-;-I))' 


( 6 . 21 ) 


( 6 . 22 ) 

(6.23) 

(6.24) 


Note that T{m + ^) = where {Q,„ := ^{^ + 1) ■ ■ ■ [^ + m - 1). By using this identity, 

we get 


(ft-l + 7(N-7))-- - (1 + 7(N-7)) 


r(ft + 7(N-7)) 
r(l + 7 (N- 7 )) 


and 


(ft - 1 + ^ + 7(2N - 7 - 1)) ■ • ■ (1 + i6 + 7(2N - j - 1)) 

It follows that 


r(ft + ^ + 7(2N-7-l)) 
r(l + /3 + 7(2N-7-l))- 


N 

SN(a,)6,7) = SN(l,j6,7)n 


M 


r(ft + 7(N - 7))r(l + |6 + 7(2N -7-1)) 
r(ft + |6 + 7(2N - 7 - l))r(l + 7(N - /)) 


. .. , A r(i + /3 + 7 ( 2 N- 7 -i)) \ N 

Nf 7 )n r(^ + L I r(ft + /i + 7 ( 2 N -y - 1 )) 


As Sj^{ci,fi, 7 ) = SN(/3,ft, 7 ) by the fact that A(1 — x) = ±A(ft), that is, S]v(a,j 6 , 7 ) is a symmetric 
function of ft and /3, at the same time, the following factor is also symmetric in ft and /3: 


^ r(ft + 7(N-7))r(/3 + 7(N-7)) 

r(ft + /3 + 7(2N-7-l)) 

It follows that the factor 


N 


Sn(1/^,7) it 

/=i 


r(l + /3 + 7 ( 2 N- 7 -l)) 
r(/l + 7(N-7))r(l + 7(N-7)) 


should be a symmetric function of ft and /I. But, however, this factor is independent of ft, therefore 
it should be also independent of /3 by the symmetry. Denote this factor by 0 ( 7 , N), we get 


SN(a,|6,7) 


/=i 

N 

c{7^N)YI 


7(N-7)) 


M 


r(ft + |6 + 7 ( 2 N- 7 -l)) 

r(ft + 7(7-i))r(/i + 7(y-i)) 

r(ft + ^ + 7(N + 7-2)) 


(6.25) 
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where 0 ( 7 , N) is independent of oc and 
To determine c{‘y,N), put a = /3 = 1; 


5^(1, 1 , 7 ) = [ |A(x)|^^dx 


N 


c(7,N)n 

/=1 


r(l + 7 (;-l))^ 
r(2 + 7 (N + y-2))- 


(6.26) 


Let y be the largest of the xi,... ,x^ and replace the other Xj by Xj = ytj, where 0 ^ fy ^ 1. 
Without loss of generality, we assume that y = Then Xj = ytj for y = 1,..., N — 1. Then 


N-l 

I A(xi,..., Xn) I'"' = • I A(fi,..., fM-i) 1 '^ ■ n (1 - 

7=1 


and the Jacobian of change of variables is 


h h ■ ■ ■ In-i 1 

1 / 0 0 ■ ■ ■ 0 

0 y 0 ■ ■ ■ 0 

0 0 • • • y 0 

Now we have 


det ( 


(6.27) 


(6.28) 


Sn( 1,1,7) = N! / |A(x)|"^dx 

r\ p ^ 1 

= Nl yT'^(^~^)y^“My- / | A(fi,..., fN_i) 1^"^ (1 - fy)^'^dfi... dlw-i 

■Jo 2o<ti<---<t]v/-isSl j^i 

“ (^i) W''~'W'dy) 0(1 “ 

that is 


SN(hl,7) = - 1 ) + 1 

c( 7.N-1) r(l + 7(; - 1))T(27 + 1 + 7(; - 1 )) 
7(N-1) + 1 r(l + 27 + l + 7(N-l+y-2)) 

c{7,N-1) ^1 t(1 + 7(;-1))T(1 + 7 + 7;) 

7(N-1) + 1 r (2 + 7(N + y-l)) 


Thus 


(6.29) 

(6.30) 

(6.31) 


N 


c( 7 .N)n 

/=i 


r(i + 7(;-i))^ 

r(2 + 7(N + 7-2)) 


Sn(1/1,7) 


^i T(l + 7(;-l))r(l + 7 + 7y) 
7(N-1) + 1 T (2 + 7(N + y-l)) 
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This implies that 


c( 7 ,N) _r(l + 7 N) 
c( 7 ,N- 1 ) “ r(l + 7 ) ' 


or 


c( 7 ,N) = c(7,l) • 
Finally we get 


c{l,N) 

c(7,N- 1) 


^(7.2) \ 

c(7.i)y 


c(7.1) ■ 


y r(i + 7 N) 
V r(i + 7) 


N 


c{-f,N) = c( 7 ,l) -n 

/=2 


r(i + 77) 
r(i + 7) 


fr r(i + 77) 

y r(i + 7)' 


where the fact that 0 ( 7 , 1 ) = 1 is trivial. 


r(i + 72) \ 
r(i + 7) y ■ 


(6.32) 


(6.33) 


(6.34) 

□ 


Remark 6.2. Note that the conclusion is derived here for integers a, /3 and complex 7 . 


A slight change of reasoning due to Askey gives (I6.25|) directly for complex a, and 7 as fol¬ 
lows. (I 6 . 88 D and the symmetry identity that Sjv(a, 7 ) = give the ratio of S]v(a,/3, 7 ) 

and SN(a, j 6 -|- m, 7 ) for any integer m, 


SN(a,j6,7) SN{a,^ + m-l,'Y) Sn{oc,^,j) 

(/3 + m-l)+ 7 (N-;) ^ i 6 + 7 (N-;) 

fJ[Cc + {^ + m-l) + 7(2N - 7 - 1) y a + ^ + 7(2N - ; - 1) 

_ rr (/I + 7 (^ ~ i))m 

“}i(« + /3 + 7 ( 2 N-;-l)),' 


(6.35) 

(6.36) 

(6.37) 


where we have used the notation 


(®)m 


T[a + m) 
T{m) ' 


m ^ 0; 


i.e. (fl)o = 1 and {a)m ■= a{a + 1) ■ ■ ■ {a + m — 1) for m ^ 1. Now 


S]v(a,j6 + m,7) 



(6.38) 

(6.39) 


Thus 


SN{a,^,y) = SN{a,^ + m,7) ■ n ^ 

' (/3 + 7(N-;■)), 


^ ^-aN-7N(N-l) 


[ 


N 

iAw7n^r‘ 




m 


;=i 


N 

n 

/=i 


{a + fi + j{2N - i - 


(/5 + 7(N-/)), 




N 

Ni 


\ ^+m-l 
X; \ 


m / 


((4:^2) 
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Denote aj = a + ^{N — j), bj = j5 + ^{N — j) and cy = a + j6 + 7(2N — j — 1), we have 

r(c+m) 

b-c _ r(&) r(c + m) _ T{b) T{m)m<= 

{b),n T{c)T{b + m) T{c)lp^- 


By using the fact that 


it follows that 


therefore 


r(m + c) 

lim —-r- 

m^oo T[m)m‘^ 


(Vc G R), 


lim 

m^co 


(c)m 

(b)m 


m 


h—c 


m 

Tic)' 



(cy)m 

{bj)m 


h ■ ~ 
mi 


N 


n 



Taking m —>• oo in (I6.40D gives rise to the following: 


SN(a,^,7) = 



exp {—Xj) dxj. 


(6.43) 


(6.44) 


(6.45) 


Furthermore, 


S]v(a, j6,7) 


/ ~ Tiaj)Tibj) \ /o“ ■ ■ • /o“ I A(^) 1^^ Uf=i x‘1 ^ exp {-Xj) dXj 

r(c;) ) nf.ir(«y) 


(6.46) 


By the symmetry of a and fi, it follows that the factor 

/o°° • • • /o" I A(^) 1^^ njli exp j-Xj) dxj 
n~ir(«y) 

is a symmetric function of a and j6. Since it is independent of /I, it is also independent of a by the 
symmetry of a and fi. Thus 


/o" • • ■ Jo" I A(^) 1^^ njli xJ ^ exp j-Xj) dxj ^ /o°° • • • J" I A(x) Ylf=i xf ^ exp j-Xj) dxj 
n;=ir(fl;) " UjLinbj) 

which is denoted by c{j,N). This indicates that 


N 


N 


10 


10 


\Hx)rYl^j ^exp(-xy)dxy = c( 7 ,N)nr(c 


A r(i + 7;) fr ry, fr + A 

ijr(i + 7)|j y- 


r(i + 7) 


(6.49) 

(6.50) 
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Remark 6.3. Now we show that 


lim I = 1 (Va G R+ U {0}). 

n^oo r(n)n“= ^ ^ 


Indeed, in order to prove this fact, we need a limit representation of the gamma function given 

oo fl+M” 

by Carl Friedrich Gauss via Euler's representation of n! = ^ : 

n\n^ 

r(z) = lim — - 

^ n^ca z{z + 1) ■ ■ ■ {Z + n) 


, r(n + a) , r(a)a(a + 1) • • ■ (a + n — 1) 
lim ^ = lim - ^—7 —^^ 

n—foo {^n — 

a(a + 1) ■ • • (a + n) n 


r(n)n‘' 


= r(a) lim 


n^oo n!n'^ n + a 

= r(«) lim ita ^ 

n—>-00 nltl^ n—>-00 n OC 


= r(a) 


r(«) 


• 1 = 1 . 


(6.51) 

(6.52) 

(6.53) 

(6.54) 


Another short and elementary proof of this fact can be derived from a resulf related to inequalities 
for Gamma function ratios IfTSlI : 


Indeed, 


+ (Vfle[0,l]). 

r(x) 

r(x + a) 


r(x) 




(6.55) 


(6.56) 


as X —)■ 00 with a fixed.(This was fhe objective of WendeTs article.) To show this, first suppose 
that a. G [0,1]. Then Eq. (16.551) gives 


/ r(x^<^ 

V X/ ^ r(x)x'* ^ 


r(; 

leasing to lim^^-oo = 1 since lim^^oo (l + ^ = lim^^oo 1 = 1. Eor a > 1, the statement 

now follows from the fact that 

r(x + x)^(x^^_l)T(x + x-l) 


T(x) 


T(x) 


We are done. 


Corollary 6.4 (Laguerre's integral). By letting xy = i/y/L and f5 = L + 1 in Selberg's integral and 
taking the limit L ^ oo, we obtain 

N N 

■ I lY — t— 'VI I I I I I I I — t— ’VII 

(6.57) 

i I I -t- 'y 1 


r r^Ar oc-i , Ar(a + 7(;-i))r(i + 7;) 

I ■■■I Gwrn*, expdx,= n ——l. 
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Corollary 6.5 (Hermite's integral). By letting Xj = yj/L and a = ^ = AL^ + 1 in Selberg's integral 
and taking the limit L ^ cx3,'we obtain 


/ + 00 |• + CO T ^ / ^\ 

•■■/ |A(^)I ^riexp 

-CX) J—oo ^ ' 


dxj = 


^27^)N/2('2A)-N(7(N-1)+1)/2J-^ 

TTr~^' 

j—\ \ “1“ t) 

Remark 6 . 6 . For an integer 7 , the last equation can also be written as a finite algebraic identity. 
Firstly, we note that 

i d \ " /■+“ 


exp iy—a^x^ — 2iflxA) x”dx = 


la dA ) i-o 


exp iy—a^x^ — 2iflxA) dx (6.59) 




_ i 


i d 


exp(-A2), 


a y 2fl d A ^ 
letting A = 0 in the above reasoning, we get 

/yjexp (-aV) x"dx = ^ exp(-A^) 

we replace a by \/a in the last equation, we get 

/I (-“5 = \/f 

thus (16.5811 therefore fakes fhe form: 

. x 7N(N-1) ^ ^ g g ^27 


A=0 


a 


TT / 1 d \ , 2 n 

exp(— x^) 

2 fldxy ’ 


i d 
1^/a dx 


exp(—x^) 


x =0 


2^/^) 


n 

li^p<q^N 


dxp dx^ 


N 


exp - E ■ 


M 


= (2a) 


(6.60) 

(6.61) 

,(6.62) 

(6.63) 


-7N(N-1)/2 ^ (n)'( 5_61) 


x=0 


(xi,...,xjv)=0 


U 7 


M 


Replacing fhe exponential by its power series expansion, one notes that the term (— xjY 
gives zero on differentiation if £ < ^N{N — l)/ 2 , and leaves a homogeneous polynomial of order 
£ — jN{N — l)/2 in the variables Xi,..., xn, if -^ > jN{N — l)/2. On setting xy = 0,; = 1,...,N, 
one sees that therefore that there is only one term, corresponding to ^ = jN{N — l)/ 2 , which 
gives a non-zero contribution. So 


n 

l^p<q^N 


( d 


N 


\dXp dXqJ 


^ -.ly / N \ ^ 




iiiV- 


At 7 ! 


(6.65) 


V7=l / M 

where £ = yN{N — l)/2. If P(x) := P(xi,...,Xn) and Q(x) := Q(xi,..., x^) are homoge¬ 
neous polynomials in x := {xi,..., Xjq) of fhe same degree, fhen a little reflection shows that 
P(d/dx)Q(x) is a constant which is also equal to Q(9/3x)P(x). Thus one can interchange the 
roles of Xj and d/dxj to get 


r — 






At 7 ! 


( 6 . 66 ) 


;=i 


where £ = jN(N — l)/2. 
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Corollary 6.7. It holds that 



where 7 is non-negative integer. 

Corollary 6 . 8 . It holds that when there is no overlap in the two sets of factors, the result is 


B{K^,K2) = 


A fl K 1 +K 2 

/ n i^-xj)^{x)dx 

i=l /^iCi+l 




+ 7(N - 0) nfii(/3 + 7(N - ;■)) 


n£t'''(« + /^ + 7(2N-;c-l)) 

where Ki,K 2 ^ 0, Ki-\- K 2 ^ N, and when there is overlap 

A A fl K1 +K2— K3 

CiK,,K2,K3) = / ••■/ n 

i=l j=Ki+l-K3 


= B{K,,K2)YI 


a + + 7 (N — k — 1) 


K + f + 1 + j{2N - k - 1)' 


(6.67) 

( 6 . 68 ) 


(6.69) 

(6.70) 

(6.71) 

(6.72) 


where Ki, K 2 , ^ 0 , Ki + K 2 — i ^3 ^ N. 


Remark 6.9. Still another integral of interest is the average value of the product of traces of the 
matrix in the circular ensembles. For example. 


Sn{p,7) 


1 ( 7 !)^ 

(27r)N (N 7 )! 


pin 

[■In 

N 

2p 


L ■ 

■1 

Y2e^^ 

n 





27 


N 


k=l 


(6.73) 


is known for 7 = 1 that Sj^{p,l) gives the number of permutations of k) in which the 

length of the longest increasing subsequence is less than or equal to N. One has in particular. 


Sn(p,1) = P'; 0 ^ p ^ N. 


It is desirable to know the integrals S]v(p, 7 ) for a general 7 . 


The following short proof of Selberg's formula is from KTH . 


Anderson's proof of Selberg's Integral. Anderson's proof depends on Dirichlet's generalization of 
the beta integral given in the following: For Re(ay) > 0, 


/■■•Xrt-vr* 


Mpo---dp„_i 


n;Lir(^y) 

r(E;L,«/)' 


(6.74) 
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where V is the set pj ^ 0/E)Lo Pj — 1- The formula is used after a change of variables. To see 
this, first consider Selberg's integral, which may be written as 

Sn = nlAn{oc,^,j)-n\ T T" ■ • • T'|(/)(0) |</>(!) |^‘^ |A^ I'" dxi ■■■ dx„, (6.75) 

Jo Jo Jo 

where 0 < Xi < X 2 < ■ ■ ■ < x„ < 1, 

^(0 = rT(^ ~ ^;) = f” ~ (pn-ii” ^ + • • • + (~1)”^0 (6.76) 

7=1 


and Afp is the discriminant of (p, so that 


Arf, — 




7=1 


n 


We now change the variables from x^, .. Xn to (po, ..., (pn-i, which are the elementary symmetric 
functions of the x/s. In fact, we have: 


= J 1 ^( 0 ) 1 “ ^|^( 1 )|^ dcpodcpi ■ ■ ■ d(p„^i, 


(6.77) 


where the integration is over all points {(po,<pi/- ■ - /(pn-i) hi which the (pj are elementary sym¬ 
metric functions of Xi,.. .,x„ with 0 < Xi < • • • < x„. Indeed, it is sufficient to prove that the 
Jacobian 

^d(pf 


det 


dx 


7 J 


— A / A . 


Observe that two columns of the Jacobian are equal when x; = Xj. Thus Yli<j{^i ~ ^j) h a factor 
of the determinant. Moreover, the Jacobian and Yli<j{^i ~ ^j) homogeneous and of the same 
degree. This proves the above Jacobian. 

We make a similar change of variables in Eq. (I6.74D . To accomplish this, set 

~ ~ (0 ^ < Cl < ■ ■ ■ < c« < 1) 

7=0 

and let 


T* = { ri(^ “ ^ 7 ) • C ;-1 < Xj < Q;/ = l,...,n\ . 


(6.78) 


. 7=1 


Next, we show that: For all (p{t) = f” — (pn-it” ^ + • • • + ( —l)”^o S TA, the following map 

{(po,(pi,...,(p„-i) ^ = iP0rPir---,Pn) e 

where (p'{t) denotes the derivative of ^(f), is a bijection and pj > 0 with Y^j^o Pi = 1- 
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Observe that 


Pi 


HQ) 

HiQ) 


[Q - Xl){Q - X2) ■ ■ ■ {Q - Xn) 

iQ - ^ 0 ) • • • {Q - Q-i){Q - Q+i) • • • iQ - Q) 


> 0 


(6.79) 


since the numerator and denominator have exactly n — j negative factors. Now let (pj{t) 
By Lagrange's interpolation formula 


t-q- 


j^O 7=0 r twl 


HQ 


(6.80) 


One can directly verify this by checking that both sides of the equation are polynomials of degree 
n and are equal at n + 1 points t = ^j,j = 0,..., n. Equate the coefficients of t" on both sides to 
get E)Lo Pj = 1- Now for a given point {po, pi,...,pn) with E-Li Pj = 1 and py > 0,; = 1, ..., n, 
define (p{t) by Eq. (16.801) . The expressions 


cpiQ = pjfjiQ) = pjiQ - ^ 0 ) • • • (Q - Q-i){Q - Q+i) • • • iQ - Q) 


and 


(piQ+i) — Pj+ifj+iiQ+i) — Pj+i{Q+i — ^0) ■ ■ ■ [Q+i — Q)iQ+i — Q+i) ■ ■ ■ {Q+i — Q) 


show that ^(^;) and (p{^j+i) have different signs and (p vanishes at some point Xj+i between Q 
and Q+i- Thus (p eV. This proves the bijection. 

We can now restate Dirichlet's formula Eq. (I6.74D as: 


fl\HQ)Q’ ^ d(po ■ ■ ■ d(p„^i 


n;L 


. 7=0 


\cp'{Ci)r-kT{‘^i) 
r (e;Li <^ 7 ) 


Indeed, 


hence 


.;-i^ ^ \HQ)r' 


J'"jvPo° Vr ^■■■PQ ^dpo---dpn-i 



(6.81) 


(6.82) 

(6.83) 


We need to verify that the Jacobian 


/ 8(po,...,p„-i) \ 

. cPn.Qj 


UW'iQ)rQ 

7=0 
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that is, dpo • • ■ dp„_i = n^o 1 9'iQ) \ ^d(po--- d^„_i or d^o • • ■ d(pn-i = Il/Lo I f'iQ) | " dpo • ■ ■ dp„_i. 


Since 


the Jacobian is 




det 


( djpor 




‘‘"'(fl),-, 


det(f') 

n”=ok'K,)|^ 

/ —l) / 


n;t„' f'Wi) 


nU^'iCj) 

1 

o 


The numerator is a Vandermonde determinant and therefore the result follows: 


J-L 


Po° ^PT Mpo---dp„_i 


„ , 'fnk'(^;Oi"dpo---dp„ 

n;Lok'(^;)r^ ^ V/=o 


--r [ nk(^;)r^ ^d^o---d^«-i- 

.frrM^ri J<p(t)evty ’ ' 


n"=ok'(^;)r^ 

The final step is to obtain the (2n — 1)-dimensional integral. Let (p{t) and 0(f) be two polynomials 
such that 


n—1 


i=l ;=1 

0 < yi < xi < y 2 < • • • < < y« < 1- 


(6.84) 

(6.85) 


The resultant of (p and O, denoted R{(p, ^), is given by 


|R(<P,0)| = 

n 


flHyj) 




!e[f5-l];;e[n] 


M 


i=l 


( 6 . 86 ) 


The absolute value of the discriminant of (p can be written as \R{(p,(p')\. That is. 


n-l 


lAfi = iAwr=n'f'( 


;=i 


The (2n — 1)-dimensional integral is 

[ 10(0) 1'^“^ |0(1) \R{cp, O) d<Po • ■ ■ d(p„_2dOo ■ • • dO„_i 


= [ i<t>(o)r-‘i<i>(i)i^‘' 


flHyj] 


7-1 


d^o • • • d^„_ 2 dOo • ■ ■ do 


n—l‘ 


(6.87) 
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Here the integration is over all (p and O defined by Eq. (I6.84D . Then we show that Selberg's 
integral An satisfies the recurrence relation: 

An {OC, fi, 7) - (rr + 7, + 7,7) • (6-88) 

In fact, integrate the (2n — 1)-dimensional integral Eq. (I6.87D with respect to d^o ■ • • ^<pn -2 and 
use 0(f) instead of cp{t) in Eq. (16.8111 to get 

7-1 




i<i.(o)r-‘10(1)1^-’ 


ri'f’fw) 

7=1 


= ^|O(0)|“ ^ ^d(po---d(pn 


r(7n) 




n^'(3/7 

7=1 


d(pQ - ■■ d(p„^2d^o ■ ■ ■ dO„_i 

2^ d^o • • • d^„_i 
dd>o • ■ ■ dO„_i 

r(7)" 


T-2 


dOo • • ■ dO„_i = ; \ A„{oc,^,j). 

r(7n) 


It remains to compute Eq. (I6.87|) in another way, set (p{t) = — Xj), and 

ocq = a,aj = 7(; = 1,.. .,n - l),a„ = /S; Xq = 0,x„ = 1 
so that Eq. (I6.87D is equal to 


f 10(0)1-' 

7(^,0) 


n—l 

ri'K 

/=i 


7-1 


dOo • ■ ■ d<t>„^id(j)o ■ ■ ■ d(pn-2 


/ ni^(^7)r^ da>o ■ ■ ■ dO„_id^o • ■ ■ d^„-2- 

h<7>/‘l’);=0 


(6.89) 


Now integrate Eq. (I6.89D with respect to dOo ■ ■ ■ dO„_i and use (p instead of (p{t) in Eq. (16.811) to 
obtain 


n—1 


r(7r)r(/3)r(7) 

r(a + j6 -6 7 (n - 1)) 


n—1 


7=1 


7-2 


k'(0)r ' dcpo ■ ■ ■ d(pn- 


Since 


k'(o)l 


n—1 
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7=1 


7=1 


7=1 7=1 
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the last integral can be written as 


r(a + /3 + 7 (n-l)) 


r(rr)r(/3)r(7)"“i 



a+7-1 


(1 - 



r(^)r(/3)r(7)”-i 




r(a + /3 + 7 (n-l)) 


Equate the two different evaluations of the {2n — 1)-dimensional integral to obtain the result. 


Finally, Selberg's formula is obtained by iterating Eq. (I6.88I) (n — 1) times. 
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